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CHAPTER  1 
Introduction 

One  of  the  many  interesting  problems  related  to  missile  guidance  and  control, 
which  has  applications  in  other  areas  as  well,  is  the  technique  known  as  gain  schedul¬ 
ing.  Frequently  theoretical  results  are  far  ahead  of  applications,  the  so  called  “theory- 
practice  gap” .  The  practice  of  gain  scheduling  seems  to  reverse  this  idea  as  the  prac¬ 
tice  is  ahead  the  theory.  The  current  state  of  the  theory  and  application  of  gain 
scheduling  is  summarized  in  the  review  papers  by  Rugh  [1]  and  by  Shamma  [2]. 

Gain  scheduling  arose  naturally  in  practice  by  engineers  who  were  versed  in  the 
classical  textbook  control  system  designs  and  were  faced  with  the  problem  of  designing 
a  controller  for  a  system  which  was  not  really  a  linear  time  invariant  system.  The 
method  converts  a  non-linear  possibly  time  var5dng  system  to  a  linear  system  by 
first  linearizing  about  a  present  operating  point.  For  example,  an  autopilot  would  be 
designed  for  given  flight  conditions  of  say,  10,000  feet  altitude,  flying  at  Mach  1.8  and 
an  angle  of  attack  of  10  degrees.  A  controller  is  designed  for  a  number  of  operating 
points  across  the  flight  envelope.  The  actual  controller  would  then  be  a  function  of 
the  usual  feedback  variables  for  the  constant  operating  condition  (linearized)  design 
and  the  operating  conditions. 

For  example,  to  control  a  non-linear  system 

we  design  a  controller  for  the  linearized  plant 

d{5x)  df{xo,uo,9)^^  ,  df{xo,UQ,9)^ 

~ir~  &  du 
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in  the  form 


6u{t)  —  k{xo,uo,  6)6x. 

Then  we  attempt  to  control  the  non-linear  plant  with  the  (non-linear)  controller 

u{t)  =  uo+  k{xo,  Uo,  9){x  -  rco). 

The  nominal  values,  Uq,  9,  Xq  are  taken  as  constant  or  time  varying  by  updating  at 
discrete  sampling  times. 

As  specified,  the  controller  is  designed  as  if  the  linearization  parameters,  Xq,  Uq, 
and  9  are  constants,  so  called  “frozen  coefiicients” .  In  order  for  this  procedure  to 
stabilize  the  system,  practitioners  of  gain  scheduling  have  developed  two  “rules-of- 
thumb”  in  the  design  of  gain  scheduled  controllers. 

1.  The  scheduling  variables  should  capture  the  plant’s  nonlinearities;  and 

2.  the  scheduling  variables  should  be  slowly  varying. 

To  illustrate  the  problem  with  the  “frozen  coefficient”  approach  consider  the  time 
varying  linear  system 

If  A(t)  were  constant,  asymptotic  stability  requires  the  eigenvalues  of  A  to  have 
negative  real  parts.  However,  a  time  varying  system  may  be  unstable  even  if  A{t)  has 
eigenvalues  with  negative  real  parts.  Stability  may  be  assured  if  ||.(4(t)||  is  sufficiently 
small,  a  requirement  which  may  not  hold  in  applications.  Another  weaker  sufficient 
condition  is  that  the  symmetric  part  of  A(t),  \{A{t)  -\-  A^(i)),  has  eigenvalues  with 
negative  real  parts  for  all  t.  This  concept  is  useful  in  extending  gain  scheduling  to 
systems  which  are  not  slowly  varying  and  is  discussed  later  in  this  report. 
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If  the  usual  “frozen  coefficient”  approach  is  taken,  one  of  the  difficulties  in  appli¬ 
cations  is  the  test  to  determine  if  the  perturbations  in  the  parameters  are  sufficiently 
small  and  the  time  variation  is  sufficiently  slow.  In  spite  of  the  lack  of  theoretical  ba¬ 
sis,  gain  scheduling  is  more  or  less  a  standard  technique  for  the  control  of  non-linear 
systems. 

Recently,  Shamma  and  Athans  [2,  3]  have  begun  to  provide  the  theoretical  foun¬ 
dation  for  gain  scheduling  and  to  quantify  these  “rules-of-thumb.”  They  have  used 
Linear  Parameter- Varying  (LPV)  systems  as  an  underlying  framework  for  nonlinear 
gain  scheduled  systems.  An  LPV  system  is  a  finite  dimensional  linear  system  whose 
dynamics  depend  continuously  on  a  time- varying  exogenous  parameter  vector  d{t). 
The  value  of  9{t)  is  not  known  a  priori  but  can  be  measured  on-line.  The  only  a  priori 
knowledge  on  9{t)  is  typically  a  knowledge  of  its  range  and  an  upper  bound  on  its  time 
rate  of  change.  This  type  of  system  can  be  used  to  describe  systems  ranging  from 
missiles  whose  exogenous  parameters  include  velocity,  altitude  and  angle-of-attack  to 
chemical  processes  whose  exogenous  parameters  might  include  temperature  and  pres¬ 
sure.  Several  other  researchers  have  considered  the  control  of  LPV  systems.  Shahruz 
and  Behtash  [4]  studied  LPV  systems  with  arbitrarily  fast  parameters.  Their  work 
presented  a  method  for  stabihzing  a  time-var3dng  linear  system  by  gain  scheduling 
without  the  slowly  varying  constraint.  They  suggested  selection  of  the  gain  matrix 
K(9),  which  is  a  function  of  the  parameter  9  such  that  the  left  hand  side  of 

A{9)  -  B{9)K{9)  +  (A{9)  -  B{9)K(9)f  =  -Q{9)  (1.1) 

is  negative  definite.  Based  on  this  equation,  they  present  an  algorithm  for  computing 
the  gain  matrix  K  (9)  and  provide  a  necessary  and  sufficient  condition  (i.e.  A22  <  0 
where  A22  is  defined  in  [4])  to  show  when  it  is  possible  to  compute  a  controller  K{9) 
for  a  fixed  9.  This  condition  must  be  satisfied  in  order  to  use  their  method  to  gain 


schedule. 
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Rugh  [5]  has  commented  that  Shahruz  and  Behtash’s  algorithm  is  too  restrictive. 
The  algorithm,  stated  in  Theorem  3.8  of  [4],  gives  necessary  and  sufficient  conditions 
not  for  stability  in  general  but  for  stability  with  V{x,t)  =  x{tfx{t)  as  a  Lyapunov 
function  for  the  system. 

For  quadratic  stability  of  a  linear  parameter- varying  system,  a  more  general  single 
quadratic  Lyapunov  function  V{x,t)  =  x{tYPx{t)  with  P  >  0  can  be  introduced 
with  corresponding  changes  to  (1.1).  Note  that  P  is  constant  for  all  6.  The  reason 
quadratic  stability  is  desirable  is  that  it  gives  a  strong  form  of  robust  stability  with 
respect  to  time  varying  parameters  [6]. 

In  this  report  we  extend  Shahruz  and  Behtash’s  work  to  include  the  more  gen¬ 
eral  condition  of  quadratic  stability.  A  new  necessary  and  sufficient  condition  for 
quadratic  stability  of  LPV  systems  will  be  derived.  This  condition  based  only  on  the 
open  loop  system  parameters  and  P  will  be  expressed  as  a  Linear  Matrix  Inequal¬ 
ity  (LMI)  which  allows  for  efficient  computation  when  checking  this  condition.  This 
condition  allows  the  computation  of  the  closed  loop  quadratic  Lyapunov  function 
without  knowledge  of  the  controller.  A  couple  of  examples  are  provided  to  demon¬ 
strate  this  new  theory  and  the  problem  of  conventional  gain  scheduling  via  “frozen 
coefficient”  controllers.  A  geometric  interpretation  of  this  new  condition  will  be  pro¬ 
vided  and  for  second  order  systems  a  method  for  graphically  representing  the  set 
of  all  valid  quadratic  Lyapunov  functions  for  the  closed  loop  system  will  be  given. 
Prom  the  P  that  is  computed,  a  stabilizing  controller  will  be  derived.  This  stabilizing 
controller  and  the  resulting  scheduling  law  is  easy  to  implement,  requiring  only  a 
few  matrices  to  be  stored  in  memory.  We  then  discuss  how  our  general  condition 
for  quadratic  stabilizability  relates  to  previous  and  concurrent  work.  Barmish  [7] 
derived  an  equivalent  A22  condition  for  uncertain  linear  systems  with  time-varying 
imcertainties,  but  the  resulting  controller  was  nonlinear  and  only  for  SISO  systems. 
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The  works  of  Khargonekar,  Zhou,  Bernussou  and  others  are  briefly  discussed.  Con¬ 
current  work  by  Becker  and  Packard  is  briefly  discussed  as  well.  In  our  discussion, 
we  show  how  our  solution  to  the  A22  condition  provides  a  solution  to  their  theorem 
of  quadratically  stabilizing  an  LPV  plant  via  output  feedback  for  which  they  did  not 
provide  a  solution. 

After  discussing  quadratic  stability  we  then  extend  our  results  to  the  even  more 
general  case  of  General  Lyapunov  Stability.  Quadratic  stability  requires  that  a  single 
Lyapunov  function  be  valid  over  the  entire  parameter  space.  This  requirement  may 
be  too  restrictive  if  it  is  known  a  priori  that  the  parameters  will  not  vary  infinetly 
fast.  If  the  parameter  variations  are  bounded  and  their  rates  of  change  are  known 
then  their  may  exist  a  time- varying  Lyapunov  function  which  is  valid  over  the  entire 
parameter  space.  We  derive  necessary  and  sufficient  conditions  for  General  Lyapunov 
Stability.  We  discuss  how  to  compute  this  time-varying  Lyapunov  function.  In  fact 
it  is  just  a  natural  extension  of  quadratic  stability.  There  is  an  extra  term  in  the 
derivation  that  is  a  function  of  P, 

Also  we  extend  our  results  of  Quadratic  Stability  to  include  £2  performance.  We 
show,  for  strictly  proper  systems,  how  to  determine  if  a  controller  can  be  constructed 
such  that  the  closed-loop  system  has  an  induced  £2  norm  less  than  a  specifled  value. 
A  missile  autopilot  example  is  provided. 

In  Chapter  3,  we  review  the  work  by  Apkarian  and  Cabinet  [8].  Their  work 
provides  an  ifoo  based  gain  scheduled  controller  for  LPV  systems.  The  parameters  are 
treated  as  uncertainties  and  are  pulled  out  into  one  large  uncertainty  block.  Collecting 
the  parameters  in  this  manner  allows  the  powerful  methods  of  If 00  theory  to  be  used. 
The  controller  is  assumed  to  have  an  LFT  dependence  on  the  parameters.  The  Small 
Gain  Theorem  and  the  Scale  Bounded  Real  Lemma  are  used  to  solve  this  problem. 
Scaling  is  introduced  to  reduce  the  conservatism  of  the  Small  Gain  Theorem  but 
the  fact  that  parameters  are  real  is  not  exploited.  This  method  provides  a  nice  way 
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to  design  a  controller  that  satisfies  an  £2  performance  criteria.  During  our  review 
we  provide  a  corrected  proof  of  their  main  theorem.  This  corrected  proof  is  due 
to  Spillman  [9]  and  Apkarian  and  Cabinet  [10].  In  addition,  we  provide  a  clear 
explanation  on  how  to  construct  a  controller.  A  discussion  of  numerical  issues  related 
to  solving  for  a  controller  is  given  as  well. 

Pulling  out  the  parameters  into  one  large  “uncertainty  block”  requires  that  the 
plant  has  a  linear  fractional  dependence  on  the  uncertainty,  i.e.  the  plant  is  a  rational 
function  of  the  parameters.  It  is  important  that  this  linear  fractional  dependence  be 
realized  in  an  efficient  manner  otherwise  the  size  of  the  matrices  get  large  resulting  in 
computational  problems.  This  type  of  realization  is  directly  related  to  the  problem 
of  N-D  realization  theory.  In  Chapter  4  we  review  the  status  of  N-D  realization 
theory.  In  fact  N-D  realization  theory  has  not  progressed  very  far  due  the  lack  of  a 
minimality  test.  We  review  some  general  methods  for  obtaining  a  realization  but  these 
do  not  provide  the  minimal  realization.  The  method  by  Belcastro  [II]  does  provide 
a  small  realization.  The  others  provide  ways  of  further  reducing  the  realization  that 
is  derived. 

Kung  et.al.  [12,  13]  do  provide  a  minimality  test  for  2-D  systems  by  computing 
greatest  common  divisors.  The  problem  that  Kung  et.al.  faced  was  that  when  the 
Createst  Common  Divisor  (CCD)  is  removed  the  system  no  longer  was  in  “state- 
space”  form  and  they  did  not  know  how  to  return  the  system  back  to  “state-space” 
form.  In  this  chapter,  we  extend  the  work  of  Rosenbrock  [14]  on  system  equivalence 
to  N-D  systems  and  use  it  to  return  a  2-D  system  back  into  “state-space”  form.  We 
provide  a  CCD  Based  algorithm  for  taking  a  nonminimal  2-D  realization,  extracting 
the  CCD  and  then  returning  it  back  into  “state-space”  form  such  that  the  resulting 
realization  is  minimal/near  minimal.  A  second  2-D  method,  the  System  Equivalent 
Based  Algorithm,  does  not  require  a  realization  but  begins  with  the  plant  description 
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P.  System  equivalent  operations  are  used  to  transform  the  realization  to  a  mini¬ 
mal/near  minimal  “state-space”  form.  This  second  method  could  easily  be  extended 
to  N-D  systems. 


CHAPTER  2 

Lyapunov  Based  Gain  Scheduling 

In  this  chapter,  we  investigate  gain  scheduling  based  on  guaranteeing  that  the 
resulting  gain  scheduled  system  satisfies  a  Lyapunov  function.  Prom  the  Lyapunov 
function  both  stability  and  performance  can  be  guaranteed  over  the  entire  parameter 
space.  In  Section  2.1,  we  investigate  gain  scheduling  from  the  standpoint  of  guaran¬ 
teeing  quadratic  stability.  A  geometric  interpretation  is  then  given.  A  discussion  is 
also  provided  on  the  connection  of  this  result  with  previous  work  in  other  related  ar¬ 
eas.  Examples  are  provided  which  show  the  practical  implementation  of  this  method 
of  gain  scheduling.  The  next  section  (Section  2.2)  extends  the  quadratic  stability  re¬ 
sult  to  general  Lyapunov  stability.  Methods  for  checking  this  condition  and  designing 
a  controller  are  discussed.  The  final  section  discusses  how  to  add  /^-performance  to 
this  new  gain  scheduled  controller  design  method.  A  missile  autopilot  design  example 
is  provided. 

For  the  purposes  of  this  chapter,  the  following  LPV  system  is  considered. 

a:(t)  A(0(t))  x{t) 

2/(t)J  [c(0(t))  £)(0(t))j  [u(t) 

where  for  all  i  >  0,  x{t)  e  Tl^,  u{t)  e  7^"",  y{t)  e  TV^,  and  the  real  matrices 
A{9)  €  B{9)  G  7^”^”“,  C(9)  G  72."»^”,  D{9)  G  axe  all  functions  of  the 

parameter  vector  9{t).  Let  this  LPV  system  be  denoted  as  Q{9,  A,  B,  C,  D)  or  more 
concisely  Qb  if  the  particular  state-space  matrices  are  clear  from  the  context. 

0(i)  =  (0i(t),...,0„,(t))GR”-  (2.2) 
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The  parameter  variations  are  not  known  a  priori  but  are  known  to  belong  to  some 
set. 


ei{t)  e  n  (2.3) 

In  addition  the  system  matrices  ^(0(t)),  B{9(t)),  C(9(t))  and  D{9{t))  must  be  bounded 
for  all  9{t)  e  With  these  definitions  for  the  LPV  system  and  its  parameters,  we 
can  now  investigate  quadratic  stabilization. 

2.1  Quadratic  Stabilization 

Quadratic  stability  can  be  formally  defined  as  follows. 

Definition  2.1  (Quadratic  Stability).  The  LPV  system  given  in  (2.1)  is  quadrat- 
ically  stable  over  the  parameter  space  if  there  exist  a  P  E  P  =  P^  >  0  such 

that  for  all  9  E  ft 

PA{9)  +  A'^{9)P<0  (2.4)  . 

This  means  that  there  exists  a  single  quadratic  Lyapunov  function,  V  {x,  t)  =x'^{t)Pxit), 
such  that  its  derivative  is  negative  definite  over  the  entire  parameter  space.  If  we  could 
design  a  gain  scheduled  controller  such  that  a  single  quadratic  Lyapunov  function  ex¬ 
ists  for  the  closed-loop  system  then  scheduling  on  slowly  varying  parameters  would 
no  longer  be  a  restriction.  In  fact  the  parameters  could  vary  infinitely  fast  and  the 
system  would  remain  stable. 

To  do  this,  we  must  first  determine  if  quadratic  stability  is  possible  and  then  design 
a  controller  that  will  achieve  it.  Shahruz  and  Behtash  in  [4]  were  close  in  designing 
such  a  controller  but  they  only  considered  the  Lyapunov  function  represented  by 
P  =  I.  Restricting  a  priori  to  this  Lyapunov  function  is  conservative  in  the  sense 
that  for  P  =  /  the  derivative  of  the  Lyapunov  function  may  not  remain  negative 
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definite  over  the  entire  parameter  space  but  there  may  exist  another  P  ^  /  which 
has  V {x,  t)  <  0. 

The  following  theorem,  which  is  a  generalization  of  the  result  in  [4]  gives  necessary 
and  sufficient  conditions  for  determining  quadratic  stabilizability  for  an  LPV  system. 
As  part  of  the  result,  a  gain  scheduled  controller  is  provided. 

Theorem  2.2  (LPV  Quadratic  Stabilizability).  The  LPV  system  defined  in  (2.1) 
is  quadratically  stable  if  and  only  if 

A22i9)  =  u^ie){A{e)w  +  wA^ie))U2ie)  < o  yoen  (2.5) 

where  W  €  7?.”^"  <  0  and  the  range  0/C/2  the  null  space  of  B(9). 

Proof.  (=^)  Suppose  the  system  is  stabilizable  then  there  exists  a  state  feedback 
K{9)  such  that  the  closed-loop  system 

x{t)  =  iA{9)  -  Bi9)K{9))x{t)  +  B{9)r{t)  (2.6) 

is  stable.  The  Lyapunov  function  corresponding  to  (2.6)  would  be 

V{x,t)  =:  x(t)'^ Px{t)  (2.7) 

Vix,9,t)  =  -x'^{t)Q{t)x{t) 

=  x'^it)  {PA{9)  -  PB{9)Ki9)  +  A'^(9)P  -  K'^{9)B'^{9)P)  x{t).  (2.8) 

Therefore 

P{A{9)  -  Bi9)K{9))  -b  (A(0)  -  B{9)K{9)fP  =  -Q{9)  (2.9) 

is  negative  definite.  An  approach  is  taken  similar  to  Shahruz’s  approach  but  makes 
no  restriction  on  P  other  than  it  has  to  be  positive  definite.  Even  though  the  closed- 
loop  Lyapunov  function  is  a  function  of  the  controller  K(9),  we  will  transform  this 
Lyapunov  function  in  two  parts.  In  one  part  the  controller  K(9)  has  complete  control 
over  and  the  other  part  the  controller  has  absolutely  no  control  over.  It  is  this 
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second  part  that  is  of  interest.  If  there  exists  a  matrix  P  >  0  that  will  make  this 
uncontrollable  part  negative  definite,  then  the  system  is  quadratically  stabilizable 
and  if  not  then  it  is  not  possible  to  make  the  system  quadratically  stable  with  an 
LPV  feedback  control. 

Proceeding  with  the  proof,  we  begin  by  expanding  (2.9).  Note  that  we  will  drop 
the  notation  which  shows  functional  dependence  on  6  to  make  the  equations  easier 
to  read. 

PA-PBK  +  A^P-K'^B'^P  =  -Q.  (2.10) 

Next  take  the  singular  value  decomposition  of  B  =  U'EV^  and  substitute  this  back 
in  for  B  yielding 

PA  -  PU'BV'^K  +  A^P  -  K'^VJfU'^P  =  -Q.  (2.11) 

Recall  that  U  is  orthogonal  and  P  is  positive  definite,  therefore  {PU)~^=  U'^P~^ 
exists.  It  will  be  convenient  to  introduce  the  following  substitution 

W  =  p-\  (2.12) 

which  makes  (PU)-^  =  V^W.  Multiplying  (2.11)  on  the  left  by  iFW  yields 

U'^A  -  +  U^WA^P  -  U^WK'^VlfU^P  =  -U'^WQ.  (2.13) 

Multiplying  on  the  right  hyWU  yields 


U'^AWU  -  T.V'^KWU  +  U'^WA^U  -  =  -U'^WQWU.  (2.14) 


Let 


A  =  U'^(AW  +  WA^)U 
L  =  V'^KWU 
Q  =  U^WQWU. 
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Then  rewriting  (2.14)  yields 


A-Y.L-  L'^YF  =  -Q. 


(2.15) 

(2.16) 
(2.17) 


(2.18) 


Now  partition  U  and  E  as 


U  =  [Ui  U2]  ,  E  = 


Si 

0 


(2.19) 


where  Ui  e  U2  e  and  E,  =  diag|aii . This 

will  partition  A  as  follows 


where 


(2.20) 


An  UfiAW  +  WA^)Ui  =  ifi  € 

ii2  :=  +  W^A^)1/2  e  (2-21) 

A22  :=  L^(ATT  +  WA^)C/2  =  A^2 


In  addition,  partition  L  as 


L{e)  =  [Ln(0)  Lum- 


(2.22) 


17 


Substituting  A(9),  S  and  L  from  (2.20),  (2.19)  and  (2.22)  into  (2.18), then  we  obtain 
An{6)  -  (Ei(0)Lii(0)  +  LUe)T.^{e))  An{e)  -  Ei(0)Li2(^) 

Since  U{9)  and  W  are  nonsingular,  the  matrix  —Q(9)  =  —U^(0)WQ{9)WU{9)  is 
negative  definite  if  and  only  if  -Qi9)  is  negative  definite  which  is  our  original  assump¬ 
tion.  Multiplying  (2.23)  from  the  left  and  right  by  [0  x'^]  and  [0  respectively, 
where  x  is  any  vector  in  72-"“”“,  we  conclude  that  A<^{&)  is  negative  definite. 

(4=)  Suppose  there  exists  a  WK  >  0  such  that  A22{9)  in  (2.23)  is  negative  definite. 
Prom  (2.16),  we  have  free  choice  of  L.  Therefore  looking  at  (2.23),  if  you  choose  the 
elements  of  Ln{9)  =  [lij{9)]  G  7?."“^”“  as 

_  ^ 1  _  choice  i  =  1, . . . n„  —  1,  j  —  i  +  1,. .  .riu  (2.24) 
On 

=  i  =  (2.25) 

where  da  >  0  are  arbitrary  real  numbers  for  alH  =  1, ...  nu  and  dy  are  the  elements 
of  A  and,  in  addition,  choose  the  matrix  Li2{9)  G  to  be 

1,2(9)  =  ^l^i9)U'[(9)(A{9)W  +  WA'^(9))U2(9),  (2.26) 

then  we  obtain  from  (2.23)  that  —Q(9)  =  diag[— dn,  ...,  —dn^n.^,  ^22(9)]  which  is 
negative  definite.  ^ 

Lemma  2.3.  If  there  exists  o  IP  >  0  such  that  A22i9)  <  0,  then  a  stabilizing  con¬ 
troller,  K(9),  can  be  chosen  as 


=  -Q(9).  (2.23) 


K(9)  =  V(9)L{9)U'^(9)P 


(2.27) 
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where  P  =  W~^  and  L{0)  -  [Ln{d)Li2{0)]  with  the  elements  of  LniO)  =  [/y(0)]  e 


•^nuXriu  QJ.Q 


aij  —ajj 

(Tii 

i  <  j 

Oit  A~dii 
2aii 

i^j 

1 

i  >  j 

(2.28) 


l^.  —  ^ ^  choice  i  =  1, . . . riu  —  1,  j  =  i  +  1, . . . (2.29) 

^ii 


.  _  hi  +  du 

^1%  -  O  5  ^  Ij  .  .  .  l^U 

2(Tii 


(2.30) 


where  da  >  0  are  arbitrary  real  numbers  for  all  i  =  1, ...n«  and  Oij  are  the 
elements  of  A.  The  matrix  L\2{0)  E  7?."“^”“”“  is 


Li2(0)  =  'E^^{9)uf{e){Aie)w  +  WA^  ie))U2{9). 


(2.31) 


Proof.  Assuming  that  W  >  0  exists  such  that  ^22(0)  <  0,  then  from  the  proof  of 


Theorem  2.2  we  see  that  the  particular  choice  of  L{9)  = 


Liii9)  Li2i9) 


from  (2.28) 


and  (2.31)  will  make  the  system  stabihzable.  Recall  (2.16) 


l{9)  =  V'^i9)K{9)WU{9). 


Solving  for  K  (9)  yields 

K{9)  =  V{9)L{9)U'^(9)P 

which  is  our  stated  controller.  D 

The  controller  K{9)  places  the  controllable  poles  of  Q{t)  which  indirectly  places 
the  poles  of  Q(t).  Performance  can  be  introduced  by  proper  choice  of  Q{t).  This 
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will  be  discussed  further  in  Section  2.3.  The  .422  condition  is  an  LMI.  LMI’s  are 
convex  and  allow  for  convex  optimization  techniques  to  be  used  to  solve  them.  Since 
(2.5)  is  a  continuous  fimction  of  time  and  must  be  satisfied  at  all  times,  this  problem 
represents  a  convex  feasibility  problem  with  an  infinite  number  of  constraints.  For 
practical  reasons  the  control  engineer  will  have  to  resort  to  gridding  the  parameter 
space  and  solve  approximations  to  the  actual  condition  (2.5).  The  disadvantage  is 
that  as  the  number  of  parameters  increase  the  number  of  constraints  corresponding 
to  each  grid  point  will  increase  exponentially.  The  limit  on  the  number  of  parameters 
will  be  determined  by  the  computational  capability  that  is  available  to  the  control 
engineer  and  the  time  frame  allowed  for  the  controller  design.  Fortunately  once  the 
.422  condition  has  been  solved,  the  scheduling  law  is  simple  to  implement  and  does 
not  require  storing  a  controller  for  each  operating  condition  (i.e.  grid  point)  and  then 
having  to  interpolate  between  controllers  as  the  system  moves  around  the  parameter 
space.  The  controller  given  in  (2.27)  only  requires  storage  of  a  few  matrices  and 
knowledge  of  the  current  operating  conditions  in  order  to  be  implemented. 

Two  examples  are  now  presented  that  illustrate  a  problem  that  can  be  encoun¬ 
tered  when  gain  scheduling  with  frozen  point  controllers  but  can  be  overcome  by  the 
Lyapunov  based  gain  scheduling  method. 

Example  2.4.  Given  the  open  loop  LPV  plant 


0  1 

0 

X  = 

X  + 

[^-(l-F.50)  -.2 

1 

y  = 

1  0 

X 

-1  <  0  <  1, 


(2.32) 
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and  a  state  feedback  matrix  K  = 


h 


k2 


with 


+  0.02 

h  =  .1, 


(2.33) 


the  closed-loop  poles  mil  be  in  left  half-plane  for  all  values  of  9  as  shown  in  Figure 
2.1.  Even  though  the  poles  are  in  the  left  half-plane  the  response  of  the  system  with 


Figure  2.1:  The  closed-loop  pole  locations  for  -1  <  0  <  1  using  the  frozen  point 
controller. 


initial  condition  2:0  =  [ii]  o.'^d  with  6  =  cos(2t)  is  unstable  as  shown  in  Figure  2.2. 


Figure  2.2:  The  response  of  the  system  with  a;o  =  [1  —  1]^  and  9  =  cos(2t)  using  the 
frozen  point  gain  scheduled  controller. 


Now  to  design  a  Lyapunov  based  controller.  First  determine  a  Lyapunov  P  that 
satisfies  the  A22  condition  by  solving  (2.5)  with  an  LMI  solver.  For  this  example  the 
system  was  gridded  over  20  evenly  spaced  points  over  the  parameter  space.  SDPSOL 
[15]  was  used  and  the  following  P  was  obtained. 

1.1250  0.3750  ,  ^ 

P  =  (2.34) 

0.3750  1.1250 

The  controller  is  designed  from  Lemma  2.3 

K{e)  =  v{e)L[e)u'^{e)p.  (2.35) 

Using  (2.28)  and  (2.31)  to  compute  L  and  with  dn  =  10,  the  matrix  L  can  be  com¬ 
puted.  Figure  2.3  shows  the  closed-loop  response  using  the  Lyapunov  based  controller 
(2.27)  with  inititial  condition  Xq  =  [ii]  and  with  6  =  cos(2t).  The  Lyapunov  based 
controller  drives  the  states  to  zero  and  does  not  go  unstable  since  the  controller  was 
designed  such  that  the  closed-loop  system,  always  satisfies  the  Lyapunov  Stability  Con¬ 
dition  under  any  parameter  trajectory. 


Figure  2.3:  The  response  of  the  system  using  the  Lyapunov  based  gain  scheduled 
controller  with  inititial  condition  xq  =  [1  —  1]^  and  with  9  =  cos(2t). 
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Example  2.5.  In  this  example  consider  the  open-loop  plant 


X  = 


— 1  +  1.5cos^(^)  1  —  1.5siii(i)  cos(i) 
1  -1 


X  + 


u, 


(2.36) 


where  t  is  our  parameter.  We  need  only  consider  0  <  i  <  27r  when  designing  con¬ 
trollers.  For  the  frozen  point  controller  let 


K  = 


2  +  1.5sin(t) 


—  sin^(t)  • 


(2.37) 


With  this  controller  the  resulting  closed-loop  system  A  matrix  is 


Ad  = 


-1  +  1.5cos^(t) 

—  1  —  1.5sin(i)  cos(i) 


1  —  1.5sm(t)  cos(i) 
— 1  +  1.5sin^(t) 


(2.38) 


which  has  closed-loop  poles  located  at  —0.25  ±  ^’0.66144  for  all  time.  This  closed- 
loop  matrix  is  Vidyasagar’s  example  of  a  time-varying  plant  which  has  constant  left 
half -plane  poles  and  yet  is  unstable.  Figure  2.4  demonstrates  this  instability  in  the 
response  of  the  system  with  inititial  condition  xo  =  [  j]. 


Figure  2.4:  The  espouse  of  the  system  using  the  frozen  point  controller  with  inititial 
condition  Xq  =  [1  1]^  . 


23 


Now  to  design  a  Lyapunov  based  gain  scheduled  controller  (2.27).  First  determine 
a  Lyapunov  P  that  satisfies  the  A22  condition.  One  such  P  is 

9.0046  8.4899 

P  =  (2.39) 

8.4899  9.0046 

The  controller  is  designed  from  Lemma  2.3 

K(e)  =  v(9)L(e)u'^(e)p.  (2.40) 

Using  (2.28)  and  (2.31)  to  compute  L  and  with  dn  =  10,  the  matrix  L  can  be  com¬ 
puted.  The  response  of  the  resulting  closed-loop  system  using  the  Lyapunov  based 
controller  with  inititial  condition  Xq  =  [J]  is  shown  in  Figure  2.5.  Under  the  same 
parameter  trajectory,  the  Lyapunov  based  gain  scheduled  design  is  stable  where  as  the 
frozen  point  design  is  not. 


Timo 


Figure  2.5:  The  response  of  the  system  using  the  Lyapunov  based  controller  with 
inititial  condition  xq  =  [1  1]^. 

□ 

To  gain  a  better  understanding  of  this  new  gain  scheduling  law,  a  geometric 
interpretation  of  the  A22  condition  is  given  in  the  next  section. 
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2.1.1  Geometric  Interpretation 

The  submatrix  A22  represents  the  part  of  the  closed-loop  Lyapunov  function  for 
which  the  controller  has  no  control.  Therefore,  if  it  is  not  negative  definite,  the  deriva¬ 
tive  of  the  Lyapunov  function  will  not  be  negative  definite  and  quadratic  stability  can 
not  be  attained.  In  this  section,  as  we  investigate  the  meaning  of  A22  <  0,  we  will  be 
able  to  derive  the  set  of  all  positive  definite  W^’s  that  satisfy  A22  <  0.  This  in  turn 
also  specifies  the  set  of  all  positive  definite  P’s  which  satisfy  the  Lyapunov  equation 
(2.9)  because  P  =  We  will  be  able  to  specify  this  set  without  specification  of 
K{9),  Note  that  A22  is  only  a  function  of  the  open  loop  system  matrices  A{d)  and 
B{e). 

The  A22  condition  (2.5)  is 

{AW  +  WA^)U2  <  0.  (2.41) 

And  since  a  quadratic  constraint  depends  only  on  the  symmetric  part  of  the  matrix, 
this  is  equivalent  to 

{U2  AWU2)  z  <0  ’iz.  (2.42) 

Let  u  =  U2Z,  i.e.  u  is  a  vector  in  the  range  of  U2.  Substituting  this  in  gives 

AW V  <0  Vi)  =  U2Z.  (2.43) 

Hence  the  quadratic  cohdition-in  (2.42)  can  be  replaced  by  another  quadratic  condi¬ 
tion  in  a  larger  space  subject  to  a  linear  side  condition.  In  order  to  get  a  geometric 
interpretation,  define  v  =  A'^v  and  substitute  this  back  into  our  equation.  Therefore 

v'^Wv  <0  Vv  —  U2Z.  (2.44) 

Prom  (2.44),  we  can  determine  its  geometrical  interpretation  for  the  second  order 
case.  Equation  (2.44)  is  the  dot  product  of  two  vectors:  v  and  Wv.  Wv  is  the 


25 


Figure  2.6:  An  example  of  what  Wv  and  v  might  look  like  plotted  on  the  ellipse 
described  by  W. 

vector,  located  at  the  point  (ul,  v2),  that  is  perpendicular  to  the  ellipse  drawn  by  W 
and  passing  through  the  point  (ul,u2).  This  is  illustrated  in  Figure  2.6.  The  vector 
V  =  A'^v.is  the  derivative  of  v  at  (ul,u2).  If  the  dot  product,  iFWv,  is  negative  then 
the  angle  between  v  and  Wv  is  greater  than  90°.  This  means  that  v  points  inside  the 
ellipse  at  the  location  (ul,  v2). 

To  summarize  the  meaning  of  (2.5)  for  the  second  order  case,  we  have  that  if  you 
can  draw  an  ellipse  described  by  W  such  that  at  the  point  of  intersection  between  the 
ellipse  and  the  line  described  by  u  =  U^z  the  resulting  v  at  this  intersection  points 
inside  the  ellipse,  then  you  can  say  that  the  system  is  stabilizable.  In  other  words, 
the  positive  definite  matrix  W  which  describes  this  ellipse  satisfies  (2.5)  and  therefore 
it  also  satisfies  the  Lyapunov  equation  (2.9)  for  some  K.  Figure  2.7  gives  an  example 
picture  of  a  valid  ellipse  for  a  given  second  order  system.  Remember  that  the  ellipses 
are  drawn  in  the  phase  plane  plot  of  and  not  A.  (This  is  why  the  variable  was 
labelled  v  instead  of  x  the  state) . 
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Figure  2.7:  An  example  of  a  valid  ellipse  for  a  given  second  order  system. 


For  the  n**'  order  system  this  means  that  if  the  velocity  vector  calculated  from 
V  =  points  inside  an  ellipsoid  W  on  the  intersection  with  v  —  U%z  (i.e.  the  range 
oiU2),  then  the  system  is  stabilizable. 

It  is  instructive  to  consider  the  set  of  all  positive  definite  matrices  which  satisfy 
the  .422  condition  for  a  fixed  0.  For  the  second  order  case  we  can  take  Wn  =  1 
without  loss  of  generality.  The  condition  >  0  implies  Wn  >  0,  Wn  >  0  aiid 
WnW22  —  Wi2  >  0.  Therefore,  we  have 

-y/w^2  <  Wi2  <  (2.45) 

By  manipulation,  the  A22  condition  in  this  case  reduces  to 


<lWi2  <  bW22  "b  O 


(2.46) 
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wza 

Figure  2.8:  An  example  of  the  set  of  all  IF  >0  which  satisfy  the  .422  condition  for  a 
given  second  order  system. 

where 

o  =  +^2itt22(-^i2 +  -^2l) +W22A22) 

b  —  (A11W21W22  +  A21M22)  (2-47) 

C  =  (A21«2i  ■422^^2i^i22) 

The  shaded  region  in  Figure  2.8  represents  a  typical  plot  of  the  set  of  PF’s  which 
satisfy  the  A22  condition  for  a  second  order  LTI  plant.  This  region  is  bounded  on 
the  left  by  the  dashed  line  corresponding  to  the  positive  definite  condition  on  IF  as 
given  in  (2.45)  and  on  the  right  by  the  solid  straight  line  corresponding  to  the  A22 
condition  as  given  in  (2.46). 

While  Figure  2.8  represents  an  LTI  system,  Figure  2.9  represents  an  LPV  sys¬ 
tem.  In  particular,  it  is  the  phase  plane  plot  of  the  LPV  plant  given  in  (2.38)  with 
cos(t  -I-  7r/18) 

B  =  .  One  hundred  grid  points  were  taken.  Each  straight  line 

—  sin(t -1-  7r/18) 

represents  a  frozen  coefficient  plant  at  the  associated  grid  point.  Looking  closely  at 
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Figure  2.9:  The  set  of  all  WF  >  0  which  satisfy  the  A22  condition  for  Vidyasagar’s 
example. 

the  inequalities  that  are  associated  with  each  straight  line  to  determine  the  valid 
Lyapunov  region  and  then  taking  the  intersection  of  all  these  regions  results  in  the 
intersection  being  the  dear  oval  shape  in  the  middle  of  the  plot.  Therefore,  this 
oval  shape  represents  the  set  of  all  valid  Lyapunov  functions  for  this  LPV  plant.  By 
choosing  a  W  inside  this  region,  the  control  designer  can  stabilize  this  plant. 

2.1.2  Relation  to  Previous  Work 

Quadratic  Stabilizability  of  Uncertain  Systems 
When  the  A23,  condition  was  first  derived,  we  did  not  realize  that  the  stability  that 
we  were  guaranteeing  was  quadratic  stability.  This  fact  did  not  become  evident  until 
we  discovered  the  concurrent  work  of  Becker  et  al.  [6]  and  Packard  et  al.  [16]  which 
investigated  quadratic  stabilization  of  LPV  systems.  These  and  related  works  will  be 
discussed  later.  A  literature  survey  was  then  conducted  on  quadratic  stabilization  of 
hnear  systems.  As  it  turns  out,  the  condition  A22  <  0  is  fundamental  in  determining 
quadratic  stabilization  of  linear  systems.  Our  A22  condition  is  an  extension  to  LPV 
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systems  of  the  work  conducted  by  Barmish  [7]  on  uncertain  linear  systems  with  time- 
varying  uncertainties.  In  particular,  he  considered  the  system 

X  =  A{r{t))x  -h  B{s(t))u  (2.48) 


e{t)  = 


r{t) 

s{t) 


(2.49) 


where  r{t)  G  represents  the  model  parameter  uncertainty  and  s{t)  G  represents 
the  input  connection  parameter  uncertainty.  The  only  difference  between  this  and 
our  LPV  definition  is  that  Barmish  does  not  make  the  assumption  that  the  parame¬ 
ters  can  be  measured  on-line.  Before  Barmish’s  work  much  of  the  work  on  satisfying 
quadratic  stabilization  centered  around  the  solution  of  so  called  matching  conditions 
[17,  18].  Barmish  enlarged  considerably  on  these  conditions  and  provided  necessary 
and  sufficient  conditions  for  quadratic  stabilizability.  In  fact  his  conditions  are  iden¬ 
tical  to  .422  <  O  .for  systems  with  time- varying  uncertainties  just  stated  in  a  different 
manner.  The  problem  is  that  the  resulting  controller  is  a  nonlinear  function  of  the 
states  and  is  only  valid  for  single  input  systems. 

Upon  Barmish’s  discovery  of  these  necessary  and  sufficient  conditions  for  quadratic 
stabilizability,  many  researchers  focused  their  efforts  on  determining  the  conditions 
needed  to  allow  for  quadratic  stabilization  of  uncertain  systems  via  linear  control. 
This  type  of  stabilization  was  referred  to  as  linear  quadratic  stabilization. 

Zhou  and  Khargonekar  [19]  pursuing  the  goal  of  linear  quadratic  stabilization 
focused  their  attention  on  systems  with  norm  bounded  time-varying  uncertainties 
entering  both  the  state  matrix.  A,  and  the  input  matrix,  B.  They  described  their 
system  as  follows. 


X  =  [A  +  AA]a;  4-  [5  -f  AB]u 


(2.50) 
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[A^  AB]  =  DF{t)E  (2.51) 

where  A,  B,  D,  E  are  known  constant  matrices  and  F{t)  G  F  C  is  the  modeling 
or  parameter  uncertainty.  The  set  F  is  assumed  to  be  Lebesgue  measurable  and  is 
defined  as 


F  :=  {F(t)  I  F'^{t)F{t)  <  /}.  (2.52) 

They  showed  that  linear  quadratic  stabilization  for  this  class  of  uncertain  systems  is 
equivalent  to  the  existence  of  a  positive  definite  matrix  solution  to  a  Riccati  equation. 
An  LTI  controller  is  given  in  terms  of  the  solution  to  this  Riccati  equation. 

Khargonekar  and  Rotea  in  [20]  show  that  quadratic  stabilizability  is  equivalent 
to  linear  quadratic  stabilizability  for  systems  with  norm  bounded  uncertainties.  In 
particular,  they  consider  the  system 


(t)  =  Ax(t)  +  Bu(t)  +  pw(t) 

(2.53) 

(t)  =  Eix(t)  +  E2u(t) 

(2.54) 

(t)  =  A(t)e(t) 

(2.55) 

where  w{t)  G  7l'‘,  e{t)  G  'RP  and  the  real  matrices  A,  B,  D,  Ei  and  E2  are  known 
and  of  appropriate  dimensions.  The  matrix  valued  function  A(t)  is  assumed  to  be 
Lebesgue  measurable  and  is  defined  as 

A(t)  G  U  :=  {[/  G  I  ||C/||  <  1}.  (2.56) 

They  show  that  quadratic  stabilization  of  this  system  is  equivalent  to  the  solution  of 
a  Riccati  equation.  Prom  this  solution  a  stabilizing  LTI  controller  can  be  constructed. 
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It  is  interesting  to  note  that  by  adding  the  assumption  that  uncertainties  (or 
parameters)  can  be  measured  on-line  then  our  A22  condition  is  necessary  and  suf¬ 
ficient  to  guarantee  linear  quadratic  stabilization.  In  fact,  for  this  case,  quadratic 
stabilization  and  linear  quadratic  stabilization  are  equivalent. 

Bernussou  et  al.  in  [21]  made  different  assumptions  on  the  uncertainty  in  order  to 
guarantee  linear  quadratic  stabilization.  He  considered  the  uncertain  linear  system 

X  =  Ax  Bu  (2.57) 

where  A  and  B  are  matrices  belonging  to  uncertainty  domains  Va  and  T>b  defined 

by 

N  N 

Va:={A\A  =  Y1  Ai  >  0,  ^  Ai  =  1}  (2.58) 

i=l  i=l 

N  N 

Vb  :={B  \B^Y1  l^j  >  Z) /"j  =  !>•  (2.59) 

j=i  j=i 

This  corresponds  to  the  case  of  linear  systems  where  the  coefficients  are  known  up  to  a 
certain  precision  defined  by  bounding  them  from  above  and  below.  This  is  commonly 
known  as  interval  plants.  In  this  case  the  uncertainties  are  not  time-varying.  The 
necessary  and  sufficient  condition  derived  is  a  form  of  A22  <  0  for  interval  plants. 

AiW  +  WAT-BjE^Jl^BT<OViVj  (2.60) 

These  conditions  are  LMI’s  and  therefore  can  be  efficiently  checked.  These  conditions 
not  only  are  for  determining  the  feasibility  of  a  constant  Lyapunov  function  (i.e. 
quadratic  stability)  but  also  for  a  constant  controller.  By  requiring  a  fixed  controller, 
the  conditions  require  additional  variables  (i.e.  R)  that  have  to  be  determined.  Now 
the  condition  given  in  (2.60)  can  be  used  to  determine  if  a  constant  controller  could 
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be  used  to  quadratically  stabilize  an  LPV  system.  If  this  condition  were  satisfied 
then  gain  scheduling  would  not  be  needed. 

Becker  and  Packard  in  [22]  and  summarized  in  [23]  extend  Bernussou’s  work  to 
a  different  class  of  uncertainties.  They  provide  necessary  and  sufficient  conditions 
for  quadratically  stabilizing  uncertain  systems  via  LTI  state-feedback  where  the  un¬ 
certainty  in  the  A  and  B  matrices  can  depend  on  both  a  correlated  uncertainty  and 
an  uncorrelated  uncertainty.  The  uncorrelated  uncertainty  is  a  special  form  of  lin¬ 
ear  fractional  uncertainty  on  real  parameters.  The  correlated  uncertainties  lie  within 
a  prescribed  convex  polytope  while  the  uncorrelated  uncertainty  lie  within  the  unit 
cube  having  dimension  equal  to  the  number  of  uncorrelated  uncertainties.  The  nec¬ 
essary  and  sufficient  conditions  are  essentially  (2.60)  checked  at  all  of  the  vertices  of 
the  uncertainty  space. 

Concurrent  Work  in  LPV  Gain  Scheduling 

Research  in  the  area  of  LPV  gain  scheduling  has  become  a  topic  of  very  intense 
interest  as  it  is  applicable  to  so  many  areas  especially  to  the  area  of  missile  autopilot 
design.  In  particular,  contemporary  with  this  research,  has  been  the  work  conducted 
by  Becker,  Packard  and  others  who  have  collaborated  with  them  in  their  research. 
Several  important  papers  in  this  area  for  which  Becker  and  Packard  were  primary 
authors  are  [6, 16,  24,  22,  23]  along  with  Becker’s  dissertation  [25]  which  is  a  summary 
for  most  of  this  work.  . 

In  the  first  part  of  [6],  Becker  et  al.  provide  necessary  and  sufficient  conditions 
for  quadratic  LPV  stabilization  via  output  feedback  with  a  dynamic  controller.  The 
LPV  plant  under  consideration  is 


x(t) 

AW  Bm 

x(t) 

y(i) 

1 

O 

_ 1 

u(t) 

(2.61) 
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The  authors  state  that  this  plant  is  quadratically  stabilizable  via  output  feedback  if 
and  only  if  the  following  two  conditions  hold: 

1.  there  exists  s.  Pk  e  Pk  =  Pk  >  0,  and  a  continuous,  bounded  function 

F  :Tl^  ^  fQj.  within  our  parameter  set, 

[A{e)  -  B{e)F{9)YPK  +  PK[A{e)  -  B[6)K{e)]  <  o,  (2.62) 

2.  there  exists  a  Pl  €  'RF^'^,Pl  =  Pf  >  0,  and  a  continuous,  bounded  function 

L  :  P*  ^  such  that  for  all  9{t)  within  our  parameter  set, 

[A{9)  -  L{9)C{9)fPL  +  Pl[A{9)  -  L{9)C{9)]  <  0.  (2.63) 


If  the  functions  F{9)  and  L{9)  exist  as  in  the  above  conditions  then 


K{9)  = 


A{9)  +  B{9)F{9)  -  L{9)C{9)  -L{9) 
F{9)  0 


(2.64) 


is  the  quadratically  stabilizing  parameter  dependent  output-feedback  controller  ex¬ 
pressed  in  compact  matrix  form.  Further  more  they  proceed  to  parameterize  all 
stabilizing  output-feedback  controllers.  What  Becker  et  al.  do  not  show  in  this  paper 
or  in  Becker’s  dissertation  [25]  is  how  to  determine  if  F{9)  and  L{9)  exist  and  if  they 
do  exist  then  how  to  compute  them.  Our  Theorem  2.2  provides  the  answer  to  both 
the  existence  and  computation  of  the  functions. 

In  the  second  part  of  [6],  the  authors  add  C2  performance  to  the  quadratic  stabi- 
hzation  problem.  The  LPV  plant  considered  is  as  follows. 


x{t) 

■  A(e) 

B,(9) 

0 

B2{9) 

x{t) 

Qlit) 

C,{9) 

0 

0 

0 

Wi{t) 

Q2{t) 

0 

0 

0 

I 

W2it) 

_y(^)_ 

_  C2{9) 

0 

I 

0 

.  . 

(2.65) 
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where  9  is  the  vector  of  parameters  as  defined  in  (2.2),  a;(t)  e  72-", 
e  Tl"”  is  the  disturbance  input,  [gf(t),  ^  is  the  performance  output, 

u{t)  €  and  y{t)  6  for  alH  >  0.  Some  assumptions  were  made  in  order  to 
get  this  form  of  an  LPV  plant  that  are  restrictive.  In  Becker’s  dissertation  [25]  these 
assumptions  are  relaxed. 

Becker  et  al.  state  that  the  open-loop  system  given  in  (2.65)  is  quadratically 
stabilizable  and  satisfies  an  £2  performance  criteria  if  and  only  if  there  exist  matrices 
Xii  >  0  e  and  Fn  >  0  G  such  that  for  all  9  in  the  parameter  space 


A{9)Yr^+YnA'^(9)-B2i9)B^{9) 

YnCf[i9)  Bi{9) 

C'i(0)Fn 

-I 

0 

Blie) 

0 

-I 

A'^{9)Xn+Xi,A{9)  -  Gj(0)C2(0) 

XnBii9)  Cf{9) 

Bf{9)Xii 

-I 

0 

C^{9) 

0 

-I 

Xn  In 

In  Vll 

(2.66) 


Given  that  these  conditions  are  feasible,  the  authors  then  provide  formulas  to  con¬ 
struct  an  LPV  controller.  The  formulas  given  are  quite  involved  and  are  rather 
awkward.  In  the  second  part  of  [16],  Packard  et  al.  improve  upon  these  controller 
formulas.  In  this  reference,  Packard  et  al.  apply  the  controller  formulas  from  [26] 
to  this  LPV  problem  and  obtain  controller  formulas  which  are  much  easier  to  use. 
We  did  discover  some  typographical  errors  when  comparing  the  formulas  presented 
in  [16]  with  those  of  [26].  The  solution  of  the  LMPs  given  in  (2.66)  require  gridding 
the  parameter  space  and  solving  an  approximation  to  (2.66).  As  long  as  the  number 
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of  parameters  is  not  too  large  then  this  gridding  of  the  space  should  not  be  too  bur¬ 
densome  computationally  but  as  the  number  of  parameters  grow  the  number  of  grid 
points  grow  exponentially. 

In  Becker’s  dissertation  he  pulls  together  the  works  from  [6], [24]  and  [16].  In  addi¬ 
tion  to  this  he  provides  necessary  and  sufficient  conditions  for  quadratic  stabilizability 
for  an  LPV  system  which  are  given  in  the  following  Lemma. 

Lemma  2.6.  [25]  The  open-loop  system  given  in  (2.1)  is  quadratically  stabilizable 
over  the  parameter  space  if  there  exists  aW  >  ^  and  a  continuous,  bounded  function 
R  :  'RP'^  xTZ  -i-  such  that 

A{e)w  +  wA^{e)  -h  B2{e)R{e)  -k  R^{e)B2{e)  <  o  (2.67) 

for  all  9  in  the  parameter  space. 

Unfortunately  Becker  does  not  indicate  how  to  compute  the  function  R{6).  One 
method  would  be  to  grid  the  space  and  have  a  different  variable  Ri  at  each  grid  point. 
This  would  not  be  practical  though  due  to  the  extremely  large  number  of  variables 
that  would  have  to  be  determined.  Fortunately  the  A22  condition  that  we  developed 
determines  stabilizability  without  having  to  compute  a  stabilizing  controller.  It  uses 
not  only  a  reduced  number  of  variables  but  an  LMI  that  is  of  reduced  dimensions. 

.  -2.2  General  Lyapunov  Stability 

Quadratic  stability  may  be  too  conservative  in  that  it  requires  a  single  quadratic 
Lyapunov  function  to  be  satisfied  over  the  entire  parameter  space.  The  implication  of 
this  single  quadratic  Lyapunov  function  is  that  the  stability  of  the  LPV  system  is  not 
dependent  on  the  time  rate  of  change  of  the  parameters.  In  fact  the  parameters  can 
change  infinitely  fast  and  the  system  will  remain  stable  as  long  as  we  have  knowledge 
of  the  current  values  of  the  parameters.  A  less  conservative  form  of  stability  would  be 
General  Lyapunov  Stability.  General  Lyapunov  Stability  does  not  hold  the  Lyapunov 
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Function  constant  as  the  system  changes  but  allows  it  to  vary  with  the  system.  In 
this  case,  however,  knowledge  of  the  maximum  time  rate  of  change  of  the  parameters 
will  be  necessary  to  insure  stability. 

|0i|<0w,  l<i<np.  (2.68) 

The  following  theorem  gives  necessary  and  sufficient  conditions  for  general  Lyapunov 
stability. 

Theorein  2.7  (General  Lyapunov  Stability).  Ths  LPV  system  defined  in  (2.1) 
is  stable  (in  the  sense  of  General  Lyapunov  Stability)  if  and  only  if  there  exists  a 
differentiable  matrix  function  P{t)  >  0  with  W{t)  =  P  ^(t)  such  that  for  all  6(t) 

A22(0,t)  =  Uj{e)iA{9)W{t)  +  W{t)A^{e)+  )U2i9)<0  (2.69) 

new  term 

where  the  range  ofU2{9)  is  the  null  space  of  B'^{9). 

Proof  (=^»)Given  our  LPV  system  as  described  in  (2.1)  with  parameters  defined  by 
(2.2)  and  having  bounds  on  there  time  rate  of  change  given  by  (2.68),  for  necessity 
we  assume  there  exists  a  controller  K  {9)  such  that  the  closed-loop  system 

x{t)  =  {A{9)  -  B{e)K{e))x{t)  +  B{e)r{t)  '  (2.70) 

is  stable.  The  general  Lyapunov  function  corresponding  to  (2.70)  would  be 

V{x,t)  =  x{tfPit)x{t)  (2.71) 

V{x,9,t)  =  -x'^  {t)Q{t)x{t) 

=  x^{t){Pit)A{9)  -  P{t)B{9)K{e)  +  A'^{9)P{t)  - 

K'^{e)B'^{9)Pit)+  ^  )x{t).  (2.72) 

new  term 
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Therefore 

P{t)A{9)  -  Pit)B{e)K{9)  +  A^{9)P{t)  -  K^{9)B^{9)P{t)  +  Pit)  =  -Qit)  (2.73) 
is  negative  definite. 

Proceeding  with  the  proof  in  the  same  manner  as  for  the  proof  of  Theorem  2.2 
but,  in  addition,  carry  along  the  “new  term”,  results  in  (2.73)  becoming 

Ai9)  -  mm  -  L'^m’^io)  =  -Qity 
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where 
iii(0)  : 

Au{0)  : 
A22(0)  : 


uf{e)  I  A{0)w{t)  +  wit)A^{e)  +  I  c^i(0)  =  Aj,{e)  e 

new  term^ 

ufi9)  ( A{e)w{t)  +  wit)A^(e)  +  ^  )  U2{e)  e 

new  term; 


ul{B)  I  A{e)wit)  +  w{t)A^ie)  +  ^  I  U2{e) 

new  term; 


Partition  L{9)  as 


m  =  [Lni9)  Li2{9)], 


(2.80) 


(2.81) 


Substituting  A{9),  E((9)  and  L{9)  from  (2.79),  (2.78)  and  (2.81)  into  (2.74),  we  obtain 


=  -Q{t).  (2.82) 


iii(^)  -  (Si(0)I'ii(^)  +  An{e\  -T.^{9)Ln{e) 

AUB)-LU0)M0)  A2i{9) 

Since  U{9)  and  W(t)  are  nonsingular,  (5(i)  <  0  if  and  only  if  (5(t)  <  0  which  is  our 
assumption.  Multiplying  (2.82)  from  the  left  and  on  the  right  by  [0  x'^]  and  [0 
respectively  where  x  is  any  vector  in  77."^”"“ ,  we  conclude  that  .422  is  negative  definite. 

(<; — )  Now  suppose  that  there  exists  a  differentiable  function  W^(t)  >  0  such  that 
A22(9)  <  0,  fi:om  (2.76),  we  have  free  choice  of  L{9).  Therefore  looking  at  (2.82),  if 
you  choose  the  elements  of  Lxi{9)  =  [/y(^)]  G  as 


kW  = 


-,  lji{9)  =  your  choice  i  =  1, . . .  n„  —  1,  7  =  i  +  1, . . .  riu 

(2.83) 


7  fo\  ®ii(^)  •  _  1  „ 

2au{9) 


(2.84) 
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where  dii  >  0  axe  arbitrary  real  numbers  for  all  z  =  1, . . .  rZu,  and  dij  are  the  elements 
of  A  and,  in  addition,  choose  the  matrix  LuiO)  6  to  be 

Li2{e)  =  12^\e)U^{Ai9)W{t)  +  Wit)Ai9f+  ^  )U2{9),  (2.85) 

new  term 

then  we  obtain  -Q{t)  =  diag[-dii , . . .  -  ,  A22 (^)]  which  is  negative  definite.  □ 

Lemma  2.8.  If  there  exists  a  differentiable  matrix  W{t)  >0  with  P{t)  =  W~^{t) 
such  that  .^22(0,^)  <  0,  then  a  stabilizing  controller,  K{9),  can  be  chosen  as 

K{9)  =  V{9)L{9)\U{9f  P{t)  (2.86) 


where  L(9)  =  [Lii{9)Li2(9)]  with  the  elements  of  Ln{9)  =  [/ij(0)]  €  are 


<Tii{e) 


i  <  j, 


hj(9)  =  { 


auW+dgiO)  •  • 

2<Ti<(fl)  *  “  J 


1  i>  j 


(2.87) 


where  dii  >  0  are  arbitrary  real  numbers  for  olZ  z  =  1, . . .  n„  and  dij  are  the  elements 
of  A.  The  matrix  Li2(9)  E  is 


Li2(0)  =  i::\9)U^{A{9)Wit)  +  W(t)A(9)'^  +  ^  )U2i9).  (2.88) 

new  term 


Proof.  Assuming  that  a  differentiable  W(t)  >  0  exists  such  that  A22(9,t)  <  0, 
then  from  the  proof  of  Theorem  2.7  we  see  that  the  particular  choice  of  L(9)  = 


Ln(9)  Li2(9) 
(2.76) 


from  (2.87)  and  (2.88)  will  make  the  system  stabilizable.  Recall 


L{9)  =  V'^(9)K(9)W(t)U(9). 
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Solving  for  K{0)  yields 

K{e)  =  v{e)L{e)u'^{e)p{t) 

which  is  our  stated  controller.  □ 

With  the  addition  of  the  “new  term”,  (2.69)  becomes  rather  difficult  to  solve  since  it 
is  a  function  of  both  W  and  W.  In  order  to  solve  (2.69)  we  will  again  have  to  result  to 
solving  an  approximation  by  girding  the  parameter  space  and  solving  for  several  single 
quadratic  Lyapunov  functions  such  that  each  Lyapunov  function  covers  a  region  of  the 
parameter  space.  The  intention  is  to  solve  for  as  few  Lyapunov  functions  as  possible. 
The  most  desirable  situation  would  be  to  find  a  single  quadratic  Lyapunov  function 
that  is  valid  over  the  entire  parameter  space.  If  this  is  the  case  then  (2.69)  reduces 
to  the  single  quadratic  Lyapunov  case  (2.5).  The  controller  K{9)  given  in  (2.86)  also 
reduces  to  the  single  quadratic  Lyapunov  controller  of  (2.27).  For  the  case  where  a 
single  quadratic  Lyapunov  function  does  not  exist,  we  will  divide  the  parameter  space 
into  several  regions  with  a  single  quadratic  Lyapunov  function  for  each  region.  The 
controller  will  then  be  constructed  as  before  except  that  the  Lyapunov  function  used 
to  construct  the  controller  will  change  depending  on  the  location  of  the  LPV  plant 
in  the  parameter  space.  There  is  a  paradigm  shift  here.  This  is  analogous  to  the 
classical  gain  scheduled  method  of  solving  for  frozen  coefficient  controllers  and  then 
interpolating  the  controllers  as  the  parameters  change,  instead,  now  we  will  solve  for 
frozen  region  Lyapunov  functions  and  interpolate  these  Lyapunov  functions  to  get  the 
current  Lyapunov  function  which  will  then  be  used  to  compute  the  controller  gains 
as  given  in  (2.86).  The  following  is  an  algorithm  to  solve  (2.69). 


Algorithm  2.9. 
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1.  Determine  the  feasibility  of  (2.69)  with  a  single  quadratic  Lyapunov  function 
(i.e.  assume  W  =  0).  If  feasible  then  implement  the  controller  as  given  by 
(2.86)  with  a  fixed  W . 

2.  If  infeasible,  divide  the  parameter  space  into  two  regions.  Then  determine  the 
feasibility  of  (2.69)  with  two  W’s.  For  example,  determine  the  feasibility  of 

U^ie)iAie)Wi  +  WiA'^iO)  +  A)U2{e)  <  O  e  Region  1 

U'^{d){A{e)W2  +  W2A^{d)  +  A)U2{e)  <  O  e  Region  2 

with  Wi,W2  >  0  where  A  is  a  positive  definite  matrix  that  allows  for  W  ^  0. 
If  feasible,  proceed  to  the  next  step.  Otherwise,  continue  dividing  the  parameter 
space,  thereby  adding  more  frozen  region  Lyapunov  functions,  until  there  is  a 
solution. 


3.  Upon  finding  a  solution,  it  becomes  necessary  to  check  (2.69).  This  check  re¬ 

quires  knowledge  of  W.  The  value  of  W  can.  be  approximated  as  follows. 


W  = 


W. 


current 


w. 


previous 


At 


(2.89) 


where  Wcurrent  is  the  current  interpolated  value  of  W  (t)  and  Wprevioua  is  the 
previous  value  ofW{t).  Unfortunately,  Wprevious  may  posses  many  different 
values.  However,  from  the  knowledge  of  the  bounds  on  9  as  given  in  ( 2. 68), 
we  can  determine  the  regions  in  the  parameter  space  that  can  reach  our  current 
location  in  one  At.  From  these  regions  we  can  compute  the  various  values  of 
W,  allowing  us  to  check  (2.69)  at  each  grid  point.  If  (2.69)  is  satisfied  at  each 
grid  point  then  the  controller  (2.86)  can  be  implemented  with  the  interpolated 
values  ofW{t). 
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Figure  2.10:  An  example  on  how  to  interpolate  the  Lyapunov  Functions. 


An  advantage  of  Algorithm  2.9  is  that  it  initially  tries  to  find  a  single  quadratic 
Lyapunov  function  and  then  if  this  is  not  possible  it  naturally  converts  to  the  General 
Lyapunov  case. 

One  method  of  interpolating  the  IF’s  is  to  consider  the  distance  from  the  current 
location  in  the  parameter  space  to  the  center  of  each  region  that  can  be  reached  from 
this  location  in  one  At.  It  is  assumed  that  the  center  of  each  region  will  have  as 
its  Lyapunov  function  the  W  that  is  associated  with  that  region.  As  an  example, 
consider  the  two  parameter  space  as  depicted  in  Figure  2.10.  The  current  W.  can  be 
interpolated  as 


W  = 


1+4 


+ 


Wi  + 


1  +  4  + 


1^2  + 


1  +  4  +  4 


Ws 


(2.90) 


with  di,  ^2  and  d,3  as  shown  in  Figure  2.10.  To  insure  that  this  IF  is  a  valid  Lyapunov 
function  we  need  to  check 


U^{9){A{9)W  +  WA^{9)  +  A)U2{9)  <  0 


(2.91) 
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where  A  is  the  same  as  the  one  used  when  computing  the  frozen  Lyapunov  regions. 
If  it  is  not  satisfied  then  solve  the  following  LMI’s 

U^{e){A{9)iW  +  AW)  +  iW  +  AW)A^{9)  +  A)U2{9)  <  0 

(2.92) 

W  + AW  >0 

with  objective  function 

Uj  =  trace(AW)  (2.93) 

The  objective  function  is  used  to  keep  the  new  W  as  close  to  the  original  interpolated 
W  as  possible.  Thereafter, 

Wnew  =  W  +  AW.  (2.94) 

If  Step  3  of  Algorithm  2.9  fails,  then  before  coming  to  the  conclusion  that  the 
system  is  uncontrollable  the  control  engineer  can  try  the  following.  Recall  that  when 
we  solved  for  the  current  regional  Wi  that  we  used  a  fixed  A  for  the.  entire  parameter 
space.  Try  solving  (2.91)  with  a  new  A  for  the  offending  locations.  This  may  mean 
setting  up  a  new  region  inside  an  existing  region.  If  a  larger  A  is  possible,  it  will 
allow  for  a  larger  W. 

The  concept  now  is  not  to  grid  the  space  with  frozen  point  controllers  but,  instead, 
to  grid  the  space  with  frozen  point  Lyapunov  functions.  These  Lyapunov  functions 
can  be  stored  in  tables  as  previously  the  controllers  were  stored  and  then  interpolated 
as  the  plant  changes  with  the  parameters.  The  controller  scheduling  law  is  given  by 
(2.86).  This  method  of  scheduling  will  guarantee  general  Lyapunov  stability. 

So  far  only  stability  has  been  investigated,  now  the  addition  of  performance  will 


be  investigated. 
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2.3  Performance 

In  this  section  the  concept  of  guaranteeing  performance  via  state-feedback  subject 
to  quadratic  stabilization  is  discussed.  Several  methods  will  be  explored  including 
miniTnizing  a  cost  function,  minimizing  the  output  energy  and  finally  guaranteeing 
£2-performance  for  a  strictly  proper  system.  A  missile  autopilot  design  example  will 
be  given  to  demonstrate  guaranteeing  £2 

2.3.1  Minimizing  a  Cost  Function 

One  method  of  introducing  performance  into  the  design  of  LPV  systems  is  to 
design  a  controller  such  that  the  resulting  closed-loop  system  minimizes  a  given  cost 
function.  The  cost  function  that  we  will  use  is  the  standard  cost  function  used  in 
control  system  design 

00 

J  —  J  +  vF  R(d)u)dt  (2.95) 

0 

where  Q  =  L'^L  >  0  and  i2(0)  >  0  for  all  0  in  the  parameter  space.  By  allowing  the 
cost  function  to  change  with  the  parameters,  allows  the  control  designer  to  set  differ¬ 
ent  performance  objectives  depending  on  where  in  the  parameter  space  the  system 
is  located.  For  example,  a  missile  autopilot  designer  for  a  ballistic  missile  will  want 
to  reduce  performance  at  apogee  when  the  dynamic  pressure  is  low  and  will  want 
increased  performance  during  the  end-game  when  the  dynamic  pressure  has  come 
back  up.  Thereom  2.10  describes  how  to  design  such  a  controller. 

Theorem  2.10.  Given  the  LPV  system  as  described  by  (2.1)  with  the  magnitude  of 
the  parameters  being  bounded  and  initial  condition  Xq,  there  exists  a  controller 

K  =  (2.96) 

that  minimizes  the  cost  function  J  —  [x^Q{9)x  +  R(6)u)dt  subject  to  quadratic 

stabilization  if  and  only  if  there  exists  a  W  =  P~^  >  0  such  that  the  following  LMI 


45 


problem  is  solvable  for  all  6  in  the  parameter  space. 

maximize  x^{W)xq 
subject  to 

Ai9)W  +  WA^ie)  -  Bi9)R-^{9)B'^{e)  WL^{9) 
L{9)W  -I 

>0 


(2.97) 


Proof.  The  idea  is  to  bound  the  cost  function  with  a  Lyapunov  function.  Suppose 
there  exists  a  quadratic  Lyapunov  function  V{t)  =  x^{t)Px{t)  such  that 

P  >  0,  ^Vix{t))  <  -{x'^Qx  +  u'^Ru)  (2.98) 

at 

then  integrating  both  sides  from  t  =  0  to  t  =  T  yields 

T 

Vix(T))-yixo)<- Jix^Qx  +  u^Ru)dt,  (2.99) 

0 

Since  V{x{T))  >  0,  we  have 

T 

V{xo)  >  j {x^ Qx  +  vF Ru)dt.  (2.100) 

0 

This  equation  shows  that  the  cost  function  can  be  bounded  from  above  by  XqPxq. 
Prom  (2.98),  we  can  derive  an  equivalent  condition  that  will  be  used  to  prove  the 
theorem.  Expanding  out  ■^V{x)  with  u  =  —Kx,  (2.98)  becomes 

x'^iA  -  BKfPx  +  x'^P{A  -  BK)x  <  -{x'^Qx  +  x^K'^RKx) 

Then  collecting  terms  it  becomes 


x'^iiA  -  BKfP  +  PiA  -BK)+Q  +  K^RK)x  <  0. 


(2.101) 
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Therefore,  if  there  exists  a  P  such  that 

F  >  0,  {A-  BKfP  +  P(A-  BK)  +  Q  +  K'^RK  (2. 102) 

then  an  upper  bound  on  the  cost  function  is  XqPxq.  Equivalently,  iiW  =  P~^,  then 
(2.102)  becomes 

VF  >  0,  W{A-  BKf  +  {A-  BK)W  +  WQW  +  WK'^RKW  <  0  (2.103) 

and  the  corresponding  upper  bound  on  the  cost  function  is  x^W^^xq. 

The  existance  of  a  controller,  K,  subject  to  quadratic  stabilization  that  minimizes 
the  given  cost  function  is  equivalent  to  (2.103)  having  a  solution  such  that  x'^Wxq  is 
maximized.  Letting  Y  =  KW,  the  QMI  in  (2.103)  becomes 

WA^  +  AIT  -  Y^B'^  -BY  +  WQW  +  Y^RY  <  0  (2.104) 

(Letting  R  =  S^S  and  by  completing  the  square) 

WA^  +  AW  +  WQW  -  BR^^B'^  +  (5F  -  {S'^y^B'^f  iSY  -  (5^)-^B^)  <  0 

(2.105) 

4=^  (Since  Y  =  KW  =  R~^B'^,  the  completed  square  term  =  0) 

WA^  +  AW  +  WQW-BR-^B'^  <0  (2.106) 

(Via  Schur  Complements) 

AW  +  WA'^  -  BR-^B'^  WiF 

<  0  (2.107) 

LW  -I 

which  completes  the  proof. 

As  an  example  to  demonstrate  Theorem  2.10,  a  controller  will  be  designed  for  the 
open-loop  system  from  Example  2.5  such  that  a  given  cost  function  is  minimized. 
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EiXample  2.11.  Recall  the  open-loop  LPV  plant  from  Example  2.5 


X  = 


— 1  +  1.5  cos^(i)  1  —  1.5  siii(i)  cos(t) 
1  -1 


x  + 


u, 


(2.108) 


with  initial  condition  xq  —  [  J  ] .  We  desire  to  design  a  controller  that  minimizes  the 
cost  function 


OO 

J  =  l{x- 
0 


1  0 
0  1 


X  +  u^(.l)u)dt. 


(2.109) 


Gridding  the  parameter  space  as  before  into  20  points  and  solving  (2.97)  for  W  = 
P~^  yields 


P  = 


9.59073  1.94165 
1.94165  0.78845 


(2.110) 


Figure  2.11  shows  the  response  of  the  system  to  the  initial  conditions.  This  controller 
quickly  drives  the  states  to  zero  yielding  a  cost  of 


J  =  4.28 


As  a  point  of  comparison  if  we  applied  this  same  cost  function  to  the  closed-loop 
system  constructed  using  the  Lyapunov  based  controller  from  Example  2. 5,  the  cost 
would  be 


J  =  11.89 


Therefore,  as  expected  the  cost  is  lower  for  our  controller  designed  with  performance 
taken  into  account. 

a 

Another  method  of  interest  in  integrating  performance  into  the  Lyapunov  based 
gain  scheduling  is  to  minimize  the  output  energy.  The  next  section  covers  this  method. 


In  this  section  we  are  interested  in  bounding  the  output  energy  of  a  system.  To 
do  this  we  extend  the  results  taken  from  Boyd  et  al.  [27]  for  LTI  systems,  where 
DyyCy  =  0,  to  the  LPV  systems  with  no  constraints  on  Dy^Cy. 

Given  a  certain  initial  state,  xo,  consider  bounding  the  output  energy  of  the 
closed-loop  system 
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then  integrating  both  sides  from  t  =  0  to  t  =  T  yields 

T 

Vix{T))  -  l/(xo)  <-J  y'^{t)yit)dt.  (2.114) 

0 

Since  V{x{T))  >  0,  we  have 

T 

y(a;o)  >  J  y'^ {t)y{t)dt.  (2.115) 

0 

This  equation  shows  that  the  output  energy  can  be  bounded  from  above  by  x^Px^. 
Prom  (2.113),  we  can  derive  an  equivalent  condition  that  will  be  used  to  solve  the  LPV 
output  energy  minimization  problem.  Expanding  out  ■^V{x)  and  inserting  (2.111) 
into  (2.113)  yields 

x^A^Px  +  x^PAcix  <  —x^C^^CciyX  (2.116) 

Then  collecting  terms 

x^iAlP  +  PAci  +  Cf^Cciy)x<0.  (2.117) 

Therefore,  if  there  exists  a  P  such  that 

P>0,  AlP  +  PA,  +  C^Ca,<0,  (2.118) 

then  an  upper  bound  on  the  output  energy  is  XqPxq.  Equivalently,  MW  =  then 
(2.118)  becomes 

W>0,  WAl  +  AciW  +  WC^^CciyW  <  0.  (2.119) 

Condition  (2.119)  is  equivalent  via  Schur  complements  to 

AaW  +  WAl  WCf^^ 

W>0, 

CclyW  -I 
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Therefore,  if  W  satisfies  (2.120)  then  an  upper  bound  on  the  output  energy  is 
x'^W-^xo. 

Now  take  this  analysis  result  and  apply  it  to  designing  an  LPV  controller,  K{6), 
for  the  LPV  system 


X  =  A{9)x  +  Bu{9)u 

y  =  Cy{e)x  +  Dyu{e)u  (2.121) 


with  u  =  K(9)x  such  that  we  minimize  the  output  energy  of  the  resulting  closed- loop 
system 


X  =  {A(9)  +  Bu{9)Ki9))x 

y  =  {C,{0)  +  D,^{e)Kmx 

The  following  theorem  describes  how  to  do  this. 


(2.122) 


Theorem  2.12.  Given  an  initial  condition  xq  for  the  LPV  system  in  (2.121)  and  that 
{A{9),B(9),Cy{9))  is  minimal;  D'^{9)Dyu{9)  is  invertible  for  all  9  in  the  parameter 
space  and  (l  -Dy„(Dj’„I)j,„)"^T>J„)  =  there  exists  a  controller 

K(e)  =  +  Die,)  (2.123) 

that  minimizes  the  output  energy  y^ydt  subject  to  Quadratic  stabilization  if  and 
only  if  there  exists  a  IT  >  0  such  that  the  following  Quadratic  Matrix  Inequality 
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(QMI)  problem  is  solvable  for  all  0  in  the  parameter  space. 

rp 

maximize  Xq  W xq 
subject  to 

+  WA^  -  BuiD^uDyur^B'^ 

-WCjDyu{D^uDyur''Bl  WCf^E’^ 
-Bu{D^^Dyu)-^Dl,CyW 

ECyW  -I 

W^>0 

Proof.  The  existance  of  a  controller,  K,ioi  the  given  system  that  minimizes  the  output 
energy  subject  to  quadratic  stabilization  is  equivalent  to  (2.120)  having  a  solution 
such  that  XqWxq  is  maximized.  Substituting  Ad,  Bd,  Cd  from  (2.122)  into  (2.120) 
and  letting  Y  =  KW  the  matrix  constraint  becomes 

AW  +  WA'^-\-BuY  +  Y^BT  {CyW  +  DyuYY 
(CyW  +  DyuY)  -I  J 

(Via  Schur  complements) 

AW  +  WA^  +  BuY  +  +  [CyW  +  Dy^YfiCyW  +  DyuY)  <  0 

(Multiplying  out  and  collecting  terms) 

AW  +  WA^  +  WC^CyW  +  Y'^D'l^DyuY 

+  {WC'^Dyu  +  Bu)Y  +  Y^{Bl  +  £>^,C,TV)  <  0  (2.127) 

(By  completing  the  square) 

AW  +  WA^  +  WC^ CyW  -  {B^  +  D^,CyWf{D';^Dyu)-\B'^  +  D'^,,CyW) 

+  {DyuY  +  Dyu{Dl,Dy,,)-\Bl  +  DI^C yW)f 

{DyuY  +  Dy,,iDlM~\Bl  +  DlCyW))  <  0  (2.128) 


(2.125) 


(2.126) 


(2.124) 

<0 
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(Since  V  =  KW  =  -K(e)  =  +  D'^^CyW)),  the  completed 

square  term  =  0) 

AW  +  WX^  +  WO^CyW  -  +  Dl,C,Wf(D';.D„)-'(Bl  +  D^C^W) 

(2.129) 

4=^  (Expanding  and  collecting  terms) 

AW  +  WA^  - 

-  BXEi^Dyyr'D^CyW  +  WCj(l  -  D,„  W'  <  0  (2.130) 

.^=^  (Via  Schur  Complements  and  letting  (I  -Dyy{D^Dyy)~^D'^)  =  since  it  is 
always  positive  semidefinite.) 

AW  +  WA'^  - 

-WC^DyyiD^Dyyr'B^y  WC^E^ 

-  By(D^Dyy)-^Dl.CyW 

ECyW  -I 

which  completes  the  proof. 

Finally  we  will  look  at  adding  C2  performance  in  the  next  section. 


2.3.3  Guaranteeing  ^-Performance 

£2-Performance  criteria  have  become  a  very  popular  way  to  specify  performance. 
As  we  begin  to  investigate  designing  Lyapunov  based  gain  scheduled  controllers  uti- 
hzing  £2"P6rformance  criteria  it  is  necessary  to  define  the  induced  £2"^orm  and  C2- 
performance  for  LPV  systems.  The  following  definition  and  >C2-performance  Lemma 
are  taken  from  [25]. 

Definition  2.13  (Induced  £2”Norm  for  an  LPV  System  [25]).  Given  a  quadrat- 
ically  stable  LPV  system  Qe  with  input  u  and  output  y,  for  zero  initial  conditions, 
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define 

||^'e||md,2  :=  sup  sup  (2.132) 

parameter  II  ^11 2 

^  space  j|it||  0 

It  G  /^2 

to  be  the  induced  C^-norm. 

Therefore,  the  induced  £2-iiorm  for  an  LPV  system  is  the  largest  input  to  output 
induced  >C2-norm  over  the  set  of  all  causal  linear  operators  described  by  the  LPV 
system.  Now  ^2-performance  can  be  quantified  as  follows. 

Lemma  2.14  (£2-Performance  Lemma  [25]).  Consider  the  strictly  proper  LPV 
system  0(9,  Ad,  Bd,  Cd,  0).  If  there  exists  a  P  >  0  e  such  that  for  all  9  in 

the  parameter  space, 

Al(9)P  +  PAd{9)  +  ^C^i9)Cd(9)  +  PBd{9)B'^(9)P  <  0,  (2.133) 

then 

1.  The  function  Ad(9)  is  quadratically  stable  over  the  parameter  space. 

2.  There  exists  a  /3  <  j  such  that  ||^5||imi,2  <  fi 

Proof.  See  [25].  □ 

Lemma  2.14  is  an  analysis  tool.  That  is  given  an  LPV  system  and  that  (2.133) 
is  satisfied,  we  can  say  that  the  system  is  quadratically  stable  and  has  a  certain  C2- 
performance  associated  with  it.  As  control  designers,  we  wish  to  design  controllers 
for  LPV  systems  such  that  we  can  guarantee  £2-P6rformance  over  the  parameter 
space.  To  guarantee  ^2-performance  for  a  strictly  proper  system  we  have  the  following 


theorem. 
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Theorem  2.15.  Consider  the  open-loop  LPV  system  There  exists 

a  state-feedback  controller  K{6),  such  that  the  resulting  closed-loop  system  Gdi^,  Ad, 
Bd,  Cd,0)  is  quadratically  stable  and  has  an  induced  £2-norm 

lle»„lli«i,2<7  (2.134) 


if  there  exists  o  PF  >  0  such  that 


u^(e)(A(e)w  +  wA^mujie)  ^u^{$)wcr’(e) 
^c{e)wu2{9)  -I 


where  the  range  ofU^ifi)  is  the  null  space  of  B (9). 


(2.135) 


Proof.  Assume  that  there  exists  >  0  such  that  (2.135)  holds.  The  strictly  proper 
open-loop  system  under  consideration  is  given  by 


X  =  A{9)x B{6)u 

y  =  C{e)x  (2.136) 

The  state-feedback  control  law  u  =  r  —  K{9)x  results  in  the  following  closed-loop 
system. 


X  =  (A(g)  -  B{9)K{e))  x  -\-  B{e)  r 

Acl  ^cl 

y  =  CiO)  x.  (2.137) 

Ca 

Then  to  insure  >C2-p6rformance  and  quadratic  stability  of  the  closed-loop  system  as 
stated  in  Lemma  2.14  we  need  to  find  a  P  >  0  and  K{6)  such  that 

{A  -  BKfP  -I-  P{A  -  BK)  -F  +  PBB'^P  <  0.  (2.138) 

Note  that  for  simplicity  the  specific  dependence  on  9  has  been  dropped  from  the 
notation.  As  was  done  previously,  substitute  in  for  B  its  singular  value  decomposition, 
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U'EV'^,  and  pre-  and  post-multiply  (2.138)  by  U'^W  and  WU  where  W  =  P-K  This 
will  yield 


A  —  YjL  —  +  Qdesired  <  0 


(2.139) 


where 


A  =  U'^iAW  +  wA^p 

L  =  V^KWU 

Qdesired  =  ^U^WCACWU  +  SE^. 


7' 


(2.140) 

(2.141) 

(2.142) 


Partition  L  =  [Ln,  L12],  U  =  [Ui  1/2],'^'  = 


El 

0 


,  where  Ui  €  72-"^"“,  U2  € 


q^nx{n-nu)  ^^d  Ei  =  diag[crii, . . .  ,  G  This  will  partition  A  as  follows 


All  Ai2 

A^  A22 


(2.143) 


where 


^11 

Ai2 

A22 


=  U'({AW  +  WA^)Ui  =  Afi  e 
=  Ul{AW  +  WA^)U2  e 
=  iq{AW  +  WA^)U2  =  AI2  e 


(2.144) 


Then  (2.139)  becomes 


Aii-(EiLii+Li;Ef) 


+  -^U^WC^CWUi  +  EiEf 


Ai2  -  E1L12  +  ^UfWa^CWU2 


+  ^U^WC^CWUi  A22  +  ^UJWC^CWU2 


<  0. 


(2.145) 
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Prom  (2. 141), we  have  free  choice  of  L.  Therefore  if  you  choose  the  elements  of 
Lii{e)  =  [lij{d)]  E  to  be 

_  ^ij  ^jjlji  “h  Qdesiredij 


lij  — 


,  Iji  =  your  choice  i  =  1, . . .  -  1,  j  =  i  +  1, . . .  n„ 

(2.146) 


,  “1"  Qdesiredij  4"  , 

la  = - ^ - ,  z  =  l,...nu 


(2.147) 


where  da  >  0  are  arbitrary  real  numbers  for  alH  =  I, ..  .Uu,  A  =  [aij],and  Q desired  = 
[Qdesiredij  {&)]■  addition,  choose  the  matrix  Li2(0)  E  7?,"“  to  be 

LniO)  =  '^i\9){Uf{e){A{e)W  +  WA'^{e))U2{9)  +  Qdesiredil-  n„,n„  +  l:  n)} 

(2.148) 

then  equation  (2.145)  becomes 


til 


-d- 


22 


-d 


Tly,  Tl\i 


U^{AW  +  WA'^)U2  +  ^U^WC^CWU2 
which  is  negative  definite  if  and  only  if 

(AW  +  WA^) U2  +  ^U^WCFCWU2  <  0. 

Prom  the  Schur  complement  this  matrix  is  negative  definite  if  and  only  if 


<0 


iq(AW  +  WA^)U^  ^u^wc'^ 
^CWU2  -I 


<  0. 


(2.149) 
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Since  we  originally  assumed  that  there  existed  a  >  0  such  that  this  matrix  was 
negative  definite  and  recalling  that  P  =  Lemma  2.14  is  satisfied  thereby  guar¬ 
anteeing  that  the  closed-loop  system  will  be  quadratically  stable  and  have  an  induced 
£2-norm  less  than  7.  □ 

Lemma  2.16.  If  there  exists  aW  >  ^  with  P  =  W~^  such  that  (2.135)  is  satisfied, 
then  a  stabilizing  controller,  K{0),  that  yields  a  closed-loop  system  with  an  induced 
C2-norm 


(2.150) 


can  be  chosen  as 


K{e)  =  v{e)L{e)u^{e)p  (2.151) 

where  L{6)  =  {Li\{d)Li2ifi)\  with  the  elements  of  elements  of  Ln(e)  =  [iii(9)]  e 

'j^nuXriu 


^ij  Qdesiredi^ 


Iji  =  your  choice  i  =  1, . . .  —  1,  j  =  i+l^..,nu 


(2.152) 


j  _  "1“  Qdesiredij  ^ii  .  _  ^ 

III  —  T  ,  I  —  i, 

2au 


,  77>jj 


(2.153) 


and  where  da  >  0  are  arbitrary  real  numbers  for  alii  =  l,...n„,  .4  =  [aij],and 
Qdesired  ~  {gdesiredijififl'  Phe  matrix  Li2(^)  ^  P'  ^  is  defined  as 


Luie)  =  sri(^){i7f {e){A[9)W  +  WA^{e))U2{Q)  +  Qd.sired{l- nu,nu  +  1:  n)}. 

(2.154) 


Proof.  Assuming  that  W  >  11  exists  such  that  (2.135)is  satisfied,  then  from  the  proof 


of  Theorem  2.15  we  see  that  the  particular  choice  of  L{9)  = 


Ln{e)  LniO) 


from 
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(2.152),  (2.153)  and  (2.154)  will  result  in  a  quadratically  stable  closed-loop  system 
having  an  induced  £2-norm<  7.  Recall  (2.141) 

L{e)  =  v'^{e)K{e)wu(d). 


Solving  for  K{6)  yields 


k{9)  =  v{e)L(e)iF{e)p 

which  is  our  stated  controller.  □ 

As  an  example  of  designing  a  controller  to  guarantee  >C2-performance  we  have  the 
following  autopilot  example. 

Example  2.17.  Consider  the  autopilot  example  taken  from  [28]  which  we  have  mod¬ 
ified.  The  modifications  to  this  example  are  that  we  allow  the  speed  and  altitude  to 
vary  instead  of  being  held  constant.  The  nonlinear  missile  dynamics  are  as  follows 


fgcos{j)Z 

wtv  ^ 

(2.155) 

f'm 

^  ~  ~r~ 

^yy 

(2.156) 

where 


d  =  reference  diameter,  .75  ft 

/  =  radian-to-degrees  conversion,  ^ 

ft 

g  =  acceleration  of  gravity,  32.2  ^ 

lyy  =  pitch  moment  of  Inertia,  182.5  slug-ft^ 

^  =  CmQSd  =  pitch  moment,  ft-lb 
Q  =  dynamic  pressure  ^ 

p  =  -2.0777E-18  alt^  -h  6.3675E-13  alP  -  6.6022E-08  alt  -I-  2.3465  E-03  = 
a  polynomial  aproximation  for  air  density 

q  =  pitch  rate, 

S  =  reference  area,  .44  ft^ 
alt  =  altitude,  10, 000  -  30, 000ft 
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V  =  velocity,  2, 000  -  4, 000^ 

Wt=  weight,  450  lb 
Z  =  CzQS  =  normal  force,  lb 
a  =  angle-of-attack,  ±20  deg 
rjz  =  normal  acceleration  in  g 

The  normal  force  and  pitch  moment  aerodynamic  coefficients  are  approximated  by 

Cz  =  (l>zaa  +  bzS  (2.157) 

Cm  —  +  bmd  (2.158) 

where 


bz  =  -0.034 

bm  =-0.206 

5  =  fin  deflection,  deg 

<f>za  =  .OOOlOSa^  -  .00945|a|  -  .170 

(j)ma  =  .0002150:2  -  .00195joi  -  .051 


The  actuators  are  modeled  as  a  second  order  system  with  transfer  function 


5is)  = 


±  1.4a;aS  ±  wl 


5c(s) 


(2.159) 


where 


6c  —  cornmanded  fin  deflection,  deg 
u)a  =  actuator  bandwidth,  150  ^ 


The  autopilot  is  required  to  control  normal  acceleration. 


Vz  = 


Wt 


(2.160) 


The  acceleration  is  typically  expressed  in  g ’s  which  is  the  reason  for  dividing  the 
normal  force  by  the  missile ’s  weight  and  not  mass.  Now  normal  acceleration  is  not  one 
of  the  states  of  the  system  but  is  a  function  of  the  fin  deflection,  6,  and  angle-of-attack, 
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a.  In  order  to  design  a  state  feedback  Lyapunov  based  controller  that  has  performance 
goals  of  controlling  acceleration,  we  need  to  have  acceleration  or  something  directly 
related  to  it  as  one  of  the  states.  To  accomplish  this  we  augment  an  integrator  to 
the  acceleration  error  which  can  then  be  fed  back  to  control  the  fin  deflection.  This 


augmentation  is  shown  in  the  block  diagram  for  the  missile  dynamic  equations  in 
Figure  2.12  as  state  ^zerrorinteg-  The  state- space  plant  description  for  the  autopilot  is 


0 

1 

0 

0 

0 

0 

0 

—1.4a; 

0 

0 

0 

0 

Clbz 

0 

C!\(f>za 

1 

0 

X  + 

0 

0 

^^26771 

0 

C24^ma 

0 

0 

0 

0 

QSBz 

L  Wt 

0 

4>zotQS 

Wt 

0 

0_ 

1 

o_ 

(2.161) 


0  0  0  0  1 

QSBz  n  4>zaQS  n  rv 
Wt  .  Wt.  ^  ^ 


(2.162) 


To  begin  the  design  of  a  Lyapunov  based  gain  scheduled  controller  with  an  C2 
performance  objective,  we  solve  (2.135)  with  C{9)  —  (0000  30)  and  7  =  1  for  W  , 
yielding 


W  = 


3.1923.B+00  -2.9917£;+02 
-2.9917E-^2  l.m2E-h05 

5.4020M2  -4,2091ie-fp0 
2.960AE-m  1.7026£-rt)2 

-1.1890E-02  -1.9795£’-01 


5.4020.B-02  2.9604J5400 

-4.2091£;-f€0  1.7026£-M)2 

2.6663£;-02  -2.4629£;-01 
-2.4629E-01  1.2872£^-f01 

-6.2319^-03  -2.4599£;-62 


-1.1890£;-02 
-1.9795^-01 
-6.2319£!-03  , 

-2.4599i;-02 
3.9060£;-03 


(2.163) 


which  can  then  be  used  to  solve  for  the  controller 


K{e)  =  Vi9)L{9)U'^{e)P 


(2.164) 


Using  (2.152),  (2.153)  and  (2.154)  and  the  free  choice  elements  of  Ln  to  be  1,  the 
matrix  L  can  be  computed.  A  unit  step  response  for  the  nominal  closed-loop  system 


Figure  2.12:  Autopilot  block  diagram. 

ft 

is  given  in  Figure  2.13.  Nominal  in  the  since  that  V  =  3000y  and  alt  =  20000 
ft  during  the  step  response.  Figure  2.14  shows  the  response  of  the  autopilot  at  the 
extremes  of  the  parameter  variations.  The  performance  varies  between  the  extremes 
because  we  held  the  Lyapunov  function  P  to  be  constant  which  forces  the  Lyapunov 
Q  to  vary  and  performance  is  related  to  this  Q.  Our  optimization  on  the  inequality 
(2.135),  pushes  the  worst  case  system  to  satisfy  the  £2  performance  criteria  almost 
exactly  while  the  other  systems  do  much  better,  as  opposed  to  all  systems  satisfying 
the  criteria  exactly. 

Figure  2.15  shows  the  unit  step  response  of  the  same  autopilot  with  both  V  and 


alt  varying  wildly. 


alt  =  20000  +  10000  cos(6t) 
V  =  3000  +  1000  cos(12t) 


(2.165) 

(2.166) 
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Figure  2.13:  Unit  step  response  of  the  autopilot  with  the  nominal  system. 

The  autopilot  remains  stable  even  under  these  unrealistic  parameter  variations.  Even 
though  the  output  acceleration  rjz  varies  quite  a  bit  during  this  four  second  interval, 
its  mean  value  is  1.0093. 


□ 


Another  method  of  LPV  gMn  scheduling  is  presented  in  the  next  Chapter.  This 


method  is  Hoo  based. 


CHAPTER  3 

Hoo  Based  Gain  Scheduling 

An  alternative  to  the  single  quadratic  Lyapunov  approach  for  designing  gain  sched¬ 
uled  controllers  for  LPV  plants  which  also  has  a  firm  theoretical  basis  has  been  pro¬ 
posed  by  Packard  and  Becker  [16],  Packard  [29],  and  Apkarian  and  Cabinet  [8].  There 
are  two  central  ideas  behind  all  three  of  these  papers.  The  first  is  to  treat  the  time- 
varying  parameters  as  uncertainties  by  pulling  them  out  into  one  large  “uncertainty” 
block.  The  resulting  LTI  portion  of  the  plant  has  a  linear  fractional  dependence  on 
this  uncertainty  block.  By  collecting  the  time-varying  portion  into  an  uncertainty 
block,  the  powerful  methods  of  H^o  control  theory  can  be  used  to  design  a  controller. 
The  second  central  idea  in  these  papers  is  in  the  assumption  that  the  controller  also 
has  a  linear  fractional  dependence  on  the  parameters.  This  controller  structure  allows 
the  controller  to  have  knowledge  of  the  parameters.  In  [16]  and  [29],  the  discrete  ver¬ 
sion  of  the  LPV  system  is  discussed  and  sufficient  conditions  for  designing  an  gain 
scheduled  controller  are  given.  Apkarian  and  Gahinet  in  [8]  discuss  both  the  discrete 
case  and  the  continuous  case  under  the  framework  of  the  scaled  bounded  real  lemma 
and  provide  sufficient  conditions  for  designing  an  gain  scheduled  controller.  The 
conditions  are  only  sufficient  due  to  the  fact  that  the  small  gain  theorem  is  used  which 
does  not  take  into  account  the  realness  of  the  parameters  as  did  the  single  quadratic 
Lyapunov  approach.  In  the  sequel,  only  the  continuous  time  case  will  be  presented. 
The  material  in  this  chapter  is  found  primarily  in  [8]  and  is  included  as  background 
for  the  present  application. 
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3.1  Hca  Gain  Scheduling  Problem 

We  begin  by  specifically  defining  the  LPV  plant  and  the  parameters.  Let  6  = 
(^1)  ,  dne)  ^  be  the  measured  parameters.  The  fact  that  the  parameters 

vary  with  time  is  assumed  but  is  dropped  from  the  notation.  The  plant  is  assumed 
to  have  a  linear  fractional  dependence  on  the  parameters.  Therefore,  the  plant  can 
be  represented  by  the  upper  LFT 


where 


q 

y 


K(P(s),e) 


(3.1) 


e  =  blockdiag(«i/,i,...  (3.2) 

w  e  is  the  disturbance  input,  q  G  is  the  controlled  output,  u  6  72.'"^  is 
the  controlled  input,  and  y  G  is  the  measured  output.  If  the  plant  can  not  be 
represented  as  an  LFT,  approximation  methods  may  be  used  to  fit  the  parameter 
dependence  to  the  LFT  model.  Equation  (3.2)  describes  the  structure  of  the  set  of 
parameters.  In  other  words,  it  tells  how  many  times  each  parameter  or  parameter 
block  is  repeated.  Realizing  a  general  LPV  plant  that  has  rational  dependence  on  the 
parameters  may  require  using  multiple  copies  of  each  parameter.  Various  methods  by 
which  to  accomplish  this  using  as  few  copies  as  possible  will  be  discussed  in  Chapter 
4.  The  upper  LFT  of  the  plant  with  the  parameters  is  depicted  in  Figure  3.1.  Writing 
the  feedback  equations  for  the  diagram  in  Figure  3.1  yields 


1 

r 

"S 

-  . 

qeis) 

Pee{s)  Pfli(s)  P02{s) 

Weis) 

q{s) 

= 

Pie{s)  Pn(s)  Pnis) 

w-is) 

y(s)_ 

P2e{s)  F2i(s)  ^22(5) 

.  . 

(3.3) 
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Figure  3.1:  A  diagram  of  the  LPV  plant  with  an  upper  linear  fractional  dependence 
on  the  parameters. 

We  =  Oqe  (3.4) 


The  controller  is  assumed  to  have  a  linear  fractional  dependence  on  0. 

u  =  Fi{K(s),e)y  (3.5) 

This  form  for  the  controller  is  not  the  most  general  one  but  allows  the  gain  scheduled 
problem  to  be  converted  into  a  standard  Ho^  problem.  Connecting  the  gain  scheduled 
controller  in  (3.5)  with  the  LPV  plant  in  (3.1)  as  depicted  in  Figure  3.2  yields 

T{P,  K,  0)  =  FiiFuiP,  0),  Fi{K,  0))  (3.6) 

With  the  plant,  parameters  and  controller  defined,  it  is  now  necessary  to  define  the 
problem.  The  precise  problem  definition  is  infiuenced  by  the  method  used  in  solving 
it.  In  the  problem  at  hand,  it  is  desired  to  design  an  Hoo  gain  scheduled  controller, 
and  in  Hqo  theory  the  small  gain  theorem  plays  a  prominent  role.  A  particular  form 
of  the  small  gain  theorem  can  be  stated  as  follows 

Theorem  3.1  (Small  Gain  Theorem).  Let  M(s)  represent  the  transfer  function 
of  a  linear  time  invariant  system  and  A(i)  be  a  stable  time-varying  perturbation. 
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Figure  3.2:  A  diagram  of  the  interconnection  of  the  LPV  plant  with  the  gain  scheduled 
controller. 

The  closed-loop  system  shown  in  Figure  3.3  is  C2-stable  for  all  A  with  ||  A||  <  ^  if 
\\M{s)\\<j. 

Proof.  See  pg  45  [30].  □ 

Essentially,  what  this  means  is  that  a  feedback  loop  constructed  of  stable  oper¬ 
ators  will  remain  stable  if  the  product  of  the  operators  is  less  than  unity.  In  terms 
of  Nyquist,  it  means  that  the  magnitude  of  the  Nyquist  plot  is  never  larger  than 
unity,  thereby  making  an  encirclement  of  the  -1  point  impossible.  This  approach  is 
somewhat  conservative  because  it  assumes  arbitrary  bounded  parameter  variations. 
Taking  into  account  parameter  variations  having  a  certain  structure  by  use  of  simi¬ 
larity  scalings  reduces  the  conservatism,  however,  constraining  the  parameters  to  be 
real  will  not  be  taken  into  account. 

Taking  the  small  gain  theorem  into  consideration,  our  problem  definition  becomes 

Find  a  control  structure  K{s)  such  that  the  LPV  controller  Fi(K{s),G)  satisfies 

i)  the  closed  loop  system  T{P,K,Q)  is  internally  stable  for  all  bounded  pa¬ 
rameter  trajectories  9{t)  such  that  7^  0^0  <  1; 
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Figure  3.3:  Typical  plant  interconnection  for  use  with  the  small  gain  theorem, 
ii)  the  induced  /^-norm  of  the  operator  T{P,  K,  ©)  satisfies 

maxJ|T(P,ii:,0)|U<7.  (3.7) 

ll®lloo<~ 

The  first  condition  i)  can  be  assumed  without  loss  of  generality.  It  is  just  a  matter  of 
redefining  the  parameters  to  satisfy  this  condition.  The  condition  ii)  is  directly  from 
the  small  gain  theorein.  We  want  to  find  a  controller  structure  K{s)  such  that  the 
feedback  loop  gain  is  less  than  unity.  Redefining  the  parameters  to  satisfy  i)  allows  for 
clear  interpretation  of  knowing  what  the  desired  7  should  be.  In  order  for  the  system 
to  remain  stable  over  the  entire  parameter  space,  we  must  have  ||r(P,  iiT,  ©)||oo  <  1, 
in  other  words,  we  must  find  a  K(s)  that  results  in  7  <  1. 

Hoo  theory  is  set  up  to  find  a  controller  K(s)  that  will  stabilize  the  plant  depicted 
in  Figure  3.1  and  not  to  find  a  controller  of  the  form  Fi(K(s),  O).  Therefore,  we 
must  find  a  way  to  put  our  problem  into  that  form.  Collecting  the  parameters  into 

©  0 

one  repeated  “uncertainty”  block,  ,  and  redrawing  Figure  3.2  as  Figure  3.4 

0  © 

0  0 

results  in  the  standard  Boo  form.  This  block  repeated  structure,  ,  will  be 

0  © 
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Figure  3.4:  A  diagram  of  the  interconnection  of  the  LPV  plant  with  the  gain  scheduled 
controller. 


denoted  0  0  0.  Figure  3.4  is  represented  mathematically  as 


r(P,A,0)  =  F„  Fi{Pa,K), 


0  0 
0  0 


where  Pa  is  the  augmented  plant  and  is  defined  in  the  following  equation. 


(3.8) 


(3.9) 


In  order  to  reduce  the  conservatism  as  mentioned  earlier,  similarity  scaling  will 
be  introduced.  The  similarity  scalings  reduce  the  conservatism  due  to  the  fact  that 
the  Hoo  norm  is  not  invariant  under  similarity  scalings.  Therefore,  in  general,  there 
exists  an  L  7^  /  such  that  the  Hoo  norm  in  (3.12)  is  minimized.  Specifically  consider 
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the  set  of  positive  definite  similarity  scalings  associated  with  the  structure  ©  in  (3.2) 

L©  =  {L  >  0 :  L©  =  ©L,  V©}  C  (3.10) 

with 

i=l 

Positive  definiteness  of  the  scalings  is  assumed  without  loss  of  generality  due  to 
the  polar  decomposition  theorem  (see  for  example  [31])  and  the  fact  that  the  Hoo 
norm  is  invariant  with  respect  to  unitary  matrices.  Positive  definiteness  insures  the 
invertibility  of  the  scaling  L  with  the  use  of  a  convex  constraint.  For  the  repeated 
block  structure,  ©  0  ©,  the  scaling  set  becomes 

{\li  L2I  1 

■^'©©0  —  s  >  0 :  Li ,  Z/3  €  Leand  L2©  =  ©L2,V©  >  .  (3.11) 

In  ^  J 

Suppose  there  is  a  scaling  matrix  L  such  that  (3.11)  holds.  .Without  loss  of 
generality  we  can  insert  blocks  as  shown  in  Figure  3.5(a).  Using  (3.11),  the  top  two 
L  and  L~^  blocks  can  be  dropped  so  that  the  diagram  can  now  be  drawn  as  shown  in 
Figure  3.5(b)  which  must  be  optimized  over  all  controllers  /ir(s)  and  scaling  matrices 
L  such  that  (3.10)  holds. 

The  following  theorem,  by  the  use  of  the  small  gain  theorem,  formally  states  a 
sufficient  condition  for  solvability  of  the  gain  scheduled  control  problem. 

Theorem  3.2  ([8]).  Consider  the  open  loop  system  Pais)  with  uncertainty  structure 
©  0  ©  along  with  its  associated  similarity  scalings  I/©©©.  If  there  exist  a  scaling 
matrix  L  E  I/©®©  and  an  LTI  controller  K(s)  such  that  the  nominal  closed-loop 
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a)  b) 

Figure  3.5:  Insertion  of  scaling  matrices  into  the  plant. 

system  Fi{Pa{s),K{s))  is  internally  stable,  and  satisfies 


L2  0 

FiiPa{s),Kis)) 

L-2  0 

0  I 

0  I 

then  Fi{K(s),&)  is  a  j-suboptimal  gain  scheduled  Hex,  controller. 


(3.12) 


Proof.  The  proof  is  a  straightforward  application  of  the  small  gain  theorem.  □ 

Equation  (3.12)  is  the  scaled  small  gain  condition.  Pictorially  (3.12)  can  be  repre¬ 
sented  as  shown  in  Figure  3.6.  The  square  roots  are  introduced  here  and  at  this  point 
are  unnecessary.  There  necessity  comes  about  from  the  Scaled  Bounded  Real  Lemma 
which  will  be  used  to  solve  (3.12) 

Theorem  3.2  is  a  particular  case  of  the  general  scaled  Hoc  problem,  which  will  be 
discussed  in  the  next  section,  but  before  proceeding  it  is  necessary  that  state-space 
realizations  for  the  plant  and  the  controller  be  defined.  The  purpose  for  this  is  that 
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Figure  3.6:  A  diagram  of  the  scaled  small  gain  condition. 


very  efficient  state-space  based  LMI  algorithms  will  be  developed  to  solve  for  the 
controller  structure  K{s).  Therefore  consider  some  state-space  realization  (ideally 
minimal  in  terms  of  the  states  and  the  parameters)  of  the  LTI  plant  P{s). 


Dee  De\  De2 

c-e 

P(s)  = 

Die  Dll  Di2 

+ 

Cl 

D2B  D21  D22 

{si  -  A)-\Be  Bi  B2) 


(3.13) 


where  the  partitioning  is  conformable  to  (3.3).  The  problem  dimensions  are  given  by 


A  6  Dee  e  Ai  e  £>22  e  (3.14) 

The  following  assumptions  are  made  concerning  the  realization  of  the  plant  P{s). 
Al)  {A,  B2,  C2)  is  stabilizable  and  detectable; 

A2)  D22  =  0. 


The  first  assumption  is  a  standard  assumption  and  is  necessary  and  sufficient  to 
allow  stabilization  of  the  plant  by  dynamic  output  feedback.  The  second  assumption 
is  made  without  loss  of  generality  but  greatly  simplifies  the  calculations.  It  should 
be  noted  that  both  the  parameter  set  0  and  the  plant  transfer  function  from  the 
disturbance  w  to  the  controlled  output  q  have  been  considered  to  be  square.  If  this 
is  not  the  case  then  by  augmenting  the  plant  with  rows  and/or  columns  of  zeros  the 
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plant  transfer  function  can  be  made  square.  This  assumption  on  the  parameters  and 
plant  transfer  functions  simplifies  the  notation  and  the  manipulations  of  the  scaling 
matrices. 

A  state-space  realization  for  the  plant  Po(s)  can  be  derived  from  (3.13)  by  aug¬ 
menting  it  with  the  parameters  due  to  the  controller. 


p.(s) 


0 

0 

0 

0 

Ir 

0 

0 

Dee 

Dei 

De2 

0 

C, 

0 

Die 

Dll 

Di2 

0 

+ 

Cl 

0 

D2e 

D21 

D22 

0 

C2 

/r 

0 

0 

0 

0_ 

_0_ 

(si  -  A)-‘(0  B,  Bi  0) 


The  state-space  realization  for  the  controller  K{s)  is 


(3.15) 


with 


K{s) 


Dku 

Dk91 


Dki6 

Dnee 


-H 


Oki 

Ck6 


(si  —  Ak)  ^(Bki  Bkb), 


(3.16) 


Ak  e 


Note  that  the  order  k  of  K(s)  is  arbitrary  at  this  point. 

3.2  General  Scaled  Problem 

In  the  last  section,  the  Hoo  gain  scheduled  problem  was  set  up  to  be  a  particular 
case  of  the  general  scaled  Hoo  problem.  By  use  of  the  scaled  bounded  real  lemma, 
the  Hoo  constraints  can  be  turned  into  LMI  constraints  which  can  now  be  efficiently 
solved  [15,  32]. 

Lemma  3.3  (Scaled  Bounded  Real  Lemma  [8]).  Consider  a  parameter  struc¬ 
ture  ©,  the  associated  scaling  set  Lq  defined  in  (3.10),  and  a  square  continuous-time 
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transfer  function  T(s)  with  realizationT{s)  =  Dd+CdisI  —  Ad)  ^Bd-  The  following 
statements  are  equivalent. 

i)  Ad  is  stable  and  there  exists  L  E  Lq  such  that 

WlHDd  +  CdisI  -  Ad)-^Bd)L-^oo  <  7-  (3.17) 


a)  There  exist  positive  definite  solutions  X  and  L  E  Lq  to  the  matrix  inequality 


AlX  +  XAd 

XBd 

BIX 

—jL 

Dl 

<0 

(3.18) 

Cd 

Dd 

-7L-1 

Proof.  See  [33]  and  references  therein. 


□ 


The  scaled  bounded  real  lemma  is  used  here  to  design  an  Hqo  gain  scheduled 
controller  as  the  Lyapunov  equation  was  used  to  design  a  gain  scheduled  controller  in 
Chapter  2.  The  scaled  bounded  real  lemma  is  an  analysis  tool.  Given  a  state-space 
realization  of  a  system,  one  can  check  whether  or  not  the  system  is  stable  and  satisfies 
a  certain  H^o  constraint.  For  our  purposes  the  state-space  system,  T{s)  will  be  the 
state-space  realization  for  the  closed-loop  system.  As  was  done  previously,  necessary 
and  sufficient  conditions  can  be  derived  in  terms  of  only  the  open-loop  system  and 
X.  Let  the  open-loop  system  be  denoted  G{s)  and  have  a  state-space  realization 


a(s) 


Dn 

Dx2 

Cx 

D21 

D22 

-h 

C2 

(si  -  A)-HBi  B2) 


(3.19) 


In  (3.19)  the  matrix  A  refers  to  the  open-loop  system  where  Ad  in  (3.17)  refers  to 
the  complete  closed-loop  system  including  the  controller. 


Theorem  3.4  (Theorem  4.1  in  [8]).  With  G{s),  ©  and  Lq  defined  as  above,  Let 
Kr  and  fCs  denote  the  bases  of  the  null  spaces  of  {Bj  DJ2  0m2xpi)  ond  {C2  D^i  OpjxpJi 
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respectively.  With  this  notation,  the  suboptimal  scaled  Hoc,  problem  is  solvable  if  and 
only  if  there  exist  pairs  of  symmetric  matrices  (R,  S)  6  72.”^”  and  (L,  J)  €  such 
that 

AR  +  RA^  RCf[  Bi 

fCl  CiR  -'yJ  Dn  <  0  (3.20) 

Bf  Dl,  -tL 

A^S  +  SA  SBi  Cf 

KAs  BfS  -jL  Dji  /Cs<0  (3.21) 

Cl  Dn  -'tJ 

>  0  (3.22) 

LeLe,JeLe,LJ  =  I.  (3.23) 

Moreover,  there  exists  suboptimal  controllers  of  order  k  if  and  only  if  (3.20)-(3.28) 
hold  for  some  quadruple  (R,  S,  L,  J)  where  R,  S  further  satisfy  the  rank  constraint 

rank{I  —  RS)  <  k.  (3.24) 

Proof.  See  Apkarian  and  Cabinet  [8].  D 

In  Theorem  3.4  the  two  matrices  R  and  S  have  replaced  the  X  matrix  as  follows 

S  N  ,  R 

E  ' 

with  (N,  M)  e  {E,  F)  e  and  {S,  R)  G  7^'*^”.  Due  to  this  substitution 

and  the  introduction  of  the  nullspaces  Kr  and  Ks,  the  LMI  (3.18)  reduces  to  the  two 


R  J 
I  S 
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LMI’s  (3.20)  and  (3.21).  Note  that  the  LMI  (3.18)  contains  the  controller  parame¬ 
ters  through  the  Ad  but  the  new  LMFs  (3.20)  and  (3.21)  involve  the  open-loop  A 
matrix.  The  relationships  given  in  (3.25)  are  used  to  find  X  and  hence  the  controller 
parameters. 

The  three  LMI  constraints  (3.20)-(3.22)  are  all  convex.  The  structure  constraints 
L,  J  Lq  are  also  convex,  but  the  constraint,  LJ  =  I,  is  strongly  non-convex. 
Developing  methods  and  algorithms  to  solve  (3.20)-(3.23)  with  the  non-convex  LJ  =  I 
constraint  remains  an  important  and  open  area  of  research. 

Theorem  3.4  is  considered  to  be  the  general  scaled  Hoc  problem  since  it  reduces 
to  the  standard  Hoo  problem  when  ©  becomes  an  arbitrary  unstructured  uncertainty 
block  (i.e.  0  €  C^’’).  In  this  case  the  scaling  set  becomes  Lq  =  {I  *  Ir'.  I  E  Tl,l  > 

0}  and  L,  J  can  be  set  to  the  identity  without  loss  of  generality.  This  will  then 
immediately  satisfy  the  LJ  =  I  constraint  leaving  the  other  convex  constraints  which 
can  efficiently  be  solved  via  various  LMI  solver  utilities  such  as  SDPSOL  [15]  or  LMI- 
LAB  [32].  This  alternative  method  of  solving  the  classical  Hoc  problem  avoids  the 
problem  with  imaginary  axis  zeros  and  the  rank  deficiencies  in  D12  and  £>21.  For  a 
detailed  discussion  see  [33]. 

3.3  Solution  of  the  Hqq  Gain  Scheduled  Problem 

Fortunately,  the  particular  structure  of  Pa{s),  which  provides  the  controller  knowl¬ 
edge  of  the  parameters,  replaces  the  strongly  non-convex  constraint  L  J  =  I  with  a 
convex  constraint.  The  following  theorem  states  the  necessary  and  sufficient  condi¬ 
tions  for  solvability  of  the  Hoc  gain  scheduled  problem. 

Theorem  3.5  (Theorem  5.1  in  [8]  and  corrected  in  [10]).  Consider  an  LPV plant 
given  by  the  LFT  interconnection  (3.1)  where  P(s)  is  a  proper  continuous-time  LTI 
plant  with  minimal  realization  (3.13),  and  0  is  the  parameter  set  given  by  (3.2). 

The  corresponding  scaling  set  for  0  is  denoted  Lq  and  defined  by  (3.10).  Finally, 
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Lz,  Jz  €  Lq, 


Lz  I 

I  Jz 


>  0 


(3.30) 


Proof.  The  proof  as  well  as  the  original  statement  of  this  theorem  as  given  in  [8] 
are  not  correct  and  are  incomplete,  but  after  personal  communication  with  Capt. 
M.S.  Spillman,  [9],  the  following  proof  is  given.  The  proof  begins  by  substituting 
our  particular  plant  structure  into  the  necessary  and  sufficient  conditions  given  in 
Theorem  3.4.  These  new  conditions  will  be  shown  to  be  equivalent  to  the  conditions 
(3.27)-(3.30).  Therefore,  the  following  substitutions  are  made.  Note  X  —^Y  reads 
replace  the  matrix  X  in  the  state-space  realization  of  G{s)  with  Y. 


0 

r  -1 

,  C2  — > 

6-2 

Cg 

0 

L  J 

Bi  — >•  (0  Bg  Bi),  B2  {B2  0) 


(3.31) 


0  0  0 

1 - 

.  0 

£>11 

0  Dgg  Dg\ 

,  Di2 

De2  0 

0  Dig  Dll 

- 1 

0 

(M 

(3.32) 


0  D2g  D21 

to 

0 

_ 1 

D21  -> 

,  D22.  — > 

7  0  0 

I 

0 

0 

1 _ 

(3.33) 


The  A  matrix  remains  unchanged  and  the  nullspaces  Mr  and  Ms  become  (Bj, 

DI2,  0)  and  (6*2,  D29,  D2].,  0),  respectively.  Upon  making  these  substitutions,  (3.20) 
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Taking  the  Schur  complement  of  (3.35)  yields 


1  -T 


PMF  +  -P 
7 


which  is  equivalent  to 


p'^Ap  + 

7 


0  Be  Bl 

1 

o 

_ 1 

0  Dee  Dei 

Jl  Jz  0 

0  Die  Dll 

1 

o 

o 

Be  Bi 
Dee  Dei 
Die  Dll 


Equation  (3.37)  is  equivalent  to 
1 


P^iP  +  -P^ 
7 


J5i 

^>10  D>ii 


Jz  0 

0  I 


0  0  0 
BJ  Dl  Die 
Bl  Dl  Dl 


P  <  0  (3.36) 


BT  DL  DL 


09 

iT 


BT  Dh  D  u 


P<0 


Jz  0 
0 


(3.37) 


:][J  3 


which  is  the  Schur  complement  of  (3.27). 

The  constraint  given  in  (3.21)  can  be  shown  to  be  equivalent  to  (3.28)  in  exactly 
the  same  manner.  The  third  constraint  given  in  (3.22)  is  already  exactly  (3.29). 
The  last  constraint  LJ  =  I  can  be  converted  to  an  LMI  constraint  by  use  of  matrix 
dilation.  Note  that  LI,  L2,  J1  and  J2  no  longer  appear  in  (3.27)-(3.29).  This 
fact  makes  them  arbitrary.  This  arbitrariness  coupled  with  matrix  dilation  converts 
LJ  =  I  into 


Ls  / 
/  Jz 


>  0 


which  is  the  last  constraint  given  in  (3.30).  Therefore,  under  the  particular  structure 
of  the  Hoo  gain  scheduled  problem  the  four  constraints  given  in  (3.20)-(3.23)  have 
been  shown  to  be  equivalent  to  (3.27)-(3.30)  which  are  all  convex.  D 


81 


These  LMI  constraints  provide  only  solvability  conditions  for  an  i/oo  gain  sched¬ 
uled  controller.  The  next  step  is  to  take  the  quadruple  (i?,  5,  L3,  J3)  and  build  the 
controller  K{s). 


3.4  Computation  of  the  Controller 

The  idea  in  constructing  a  controller  is  to  first  construct  an  X  from  R  and  S  and 
an  L,  which  in  turn  also  specifies  J,  from  L3  and  J3.  Once  these  are  constructed, 
the  scaled  bounded  real  lemma  can  be  used  to  solve  for  the  controller,  which  is  now 
an  LMI  in  terms  of  the  controller.  Apkarian  and  Gahinet  [8]  showed  that  all  of  the 
controller  parameters  can  be  gathered  into  one  matrix,  12  defined  in  (3.41),  and  the 
scaled  bounded  real  lemma  could  be  rewritten  as 


X  0 
0  I 


V'^QQ  + 


X 

0 


(3.39) 


where  Q,  and  Q  are  given  by 


A^X  +  XAo 

XBo 

BIX 

-jC 

,  V  :=  (e’’,  0, 

(3.40) 

Co 

Vn 

iJ 

Ak 

Bki 

Bk9 

Cki 

Dku 

Bk19 

g  ^(fc+m2+r)x(fe+p2+r)^  q  ^  q) 

(3.4i) 

Ck9 

Dk91 

Bk99 

and  where 


An  — 


A  0 

0  Ofcxfc 


Bo  = 


0  Bq  Bi 

Ofcxr  0  0 


0  B2  0 

Ik  0  Ofexr 


0  Or 


=  Ce  0 
Cl  0 


0  0 


0  Dee  Dei 

0  Die  Dll 

0  h 
C2  0  , 

0  OrxA: 


Orxfc  0  If 

—  0  De2  0 

0  Dio  0 


0  0 


0  D25  D21  ) 


400 


L  0 


J  =  C-\ 


(3.42) 


In  (3.39)  the  unknowns  are  and  certain  elements  of  X,  C  and  J ,  the  following 
algorithm  shows  how  (3.39)  may  be  used  to  construct  a  controller. 

First  we  consider  finding  a  solution  X  from  (3.25)  given  S  and  R.  Prom  (3.25) 


MiV’’  =  I-RS. 


(3.43) 


Using  the  singular  value  decomposition  of  /  —  RS, 


I-RS  =  Mm. 

M 

By  identification  of  terms,  we  have  matrices  M  and  m  and 


(3.44) 


k  =  rank(/  —  RS). 


(3.45) 
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Again  from  (3.25)  since  M  is  full  rank,  we  have 

E  =  (3.46) 

which  completely  determines  X. 

We  also  need  L  which  may  be  found  from 

(311) 

^  Li  L2 
Ll  L3 

Ji  J2 

J  := 

Jz 

It  follows  that 

I  -  =  LIJ2  (3.48) 

and  J2  and  L2  may  be  found  from  the  singular  value  decomposition  of  I-  JzLz.  Note 
that  in  general  J3  is  a  block  diagonal  matrix  where  the  number  and  size  of  blocks 
is  determined  by  the  number  of  parameters  and  how  many  times  each  parameter  is 
copied.  Therefore,  compute  the  singular  value  decomposition  of  each  diagonal  block 
and  then  take  these  decompositions  and  build  L2  and  J2.  For  example,  let  UiEiV^ 
be  the  singular  value  decomposition  of  {I  —  Lz^Jz^)  which  corresponds  to  the  block 
of  (/  —  L3J3),  then  L2i  and  72,  can  be  computed  as 

^S^=svd(/-L3,J3j. 

Therefore  L2  =  diag(L2i , . . . , L2„^ )  and  J2  =  diag( J21  >  •  •  •  >  J2ng)-  Then  since 

L1J2  =  — T2J3, 

L,  =  -L2J^4(J,4)-' 


Lz  and  Jz  in  a  similar  manner.  Prom 


(3.47) 


=  L-‘. 


(3.49) 
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which  completes  the  determination  of  L. 

The  controller  parameters  are  the  only  remaining  unknowns  in  (3.39)  so  that 
(3.39)  may  be  solved  as  an  LMI  for  fi. 

There  are  various  other  explicit  methods  of  determining  the  controllers.  These 
methods  are  described  in  [34,  33,  35].  In  particular  in  [35],  efficient  and  numerically 
reliable  formulas  are  presented  for  both  the  full-order  and  reduced-order  cases.  These 
controllers  axe  simple  extensions  to  the  standard  Hoo  central  controllers.  The  advan¬ 
tage  of  these  controllers  are  in  their  efficiency  and  numerical  reliability.  If  there  are  no 
other  constraints  on  the  controller  or  the  closed-loop  properties  of  the  system,  then 
these  explicit  formulas  can  be  used.  If  on  the  other  hand  further  constraints  that 
can  be  expressed  as  LMI’s  are  placed  on  the  controller,  then  the  above  method  to 
compute  the  controller  should  be  used.  Computing  the  controller  as  the  solution  to 
an  LMI,  offers  the  flexibility  of  adding  more  constraints  to  the  controller.  The  disad¬ 
vantage  is  that  it  is  much  more  computationally  intensive  and  tends  to  be  numerically 
ill-conditioned.  Section  3.6  addresses  methods  of  reducing  the  numerical  problems 
related  to  the  above  method  for  solving  for  the  controller.  There  is  an  assumption 
that  the  LFT  realization  of  the  controller  always  exists  but  in  fact  it  may  not.  Both 
Apkaxian  and  Cabinet  [8]  and  Packard  [36]  discuss  this  situation  and  provide  work 
arounds  to  this  problem. 

3.5  Comparison  of  H no  and  LPV  Controllers 

It  is  important  to  also  know  whether  or  not  an  Hoo  gain  scheduled  controller 
provides  better  performance  than  a  single  LTI  robust  controller.  If  it  does  not,  then 
there  is  no  reason  to  go  to  all  the  extra  trouble  to  gain  schedule.  It  is  not  clear  that  the 
Hoo  gain  scheduled  controller  will  take  into  account  the  knowledge  of  the  parameters, 
since  Hoo  theory  is  predicated  on  reducing  the  affects  of  system  uncertainty  and  the 
system  parameters  are  treated  as  uncertainties. 
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To  see  how  the  Hqo  gain  scheduling  problem  does  in  fact  provide  better  perfor¬ 
mance,  first  consider  the  LTI  robust  control  problem.  Given  the  plant  realization 
P(s)  in  (3.13),  find  a  parameter  independent  LTI  controller  K(s)  such  that 


max  ||F«(Fi(P(s),ii:(s)),0)||oo  <  7  (3.50) 

lieiloo<i 

By  applying  the  small  gain  theorem,  a  sufficient  condition  for  solvability  is  the  exis¬ 
tence  of  a  scaling  matrix  L  E  Lq  and  of  an  LTI  controller  K (s)  such  that 


0 


0 

I 


Fi{P{s),K{s)) 


L-2  0 

0  I 


<7 


(3.51) 


Recall  that  Theorem  3.4  addressed  this  problem  and  provides  the  following  solvability 
conditions. 


Theorem  3,6.  Assuming  Al)-A2),  problem  (3.51)  is  solvable  if  and  only  if  there 
exists  a  quadruple  {R,  S,  Lz,  Jz)  satisfying  (3.27)-(3.30)  together  with 

LzJz  =  It  (3.52) 


Proof.  Same  as  the  proof  of  Theorem  3.4  except  replace  the  scalings  L  and  J  with 


Lz 

0 


‘Pi 


Jz  0 


0  J. 


Pi 


,  respectively. 


□ 


Comparing  Theorems  3.5  and  3.6,  the  only  difference  is  in  the  condition  placed  on 
Lz  and  Jz.  Theorem  3.5  states  that  Lz  and  Jz  must  satisfy  (3.30)  which  is  equivalent 
to 


Lz  ^0,  J3  >  0,  ^min{,LzJz^  ^  1. 


(3.53) 
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The  interpretation  of  this  is  that  the  set  of  admissible  scaling  matrices  is  any  positive 
scalings  L3,  J3  such  that 

XiiLzJz)>l,i  =  lr-  ,r,  (3.54) 

while  the  set  of  admissible  scaling  matrices  for  robust  LTI  controller  synthesis  is  any 
positive  scaling  L3,  J3  such  that 

Xi{LzJz)  =  l,i  =  l,---  ,r.  (3.55) 

Clearly,  the  robust  LTI  set  described  by  (3.55)  is  smaller  than  the  set  described  by 
(3.54).  In  fact,  (3.55)  is  the  boundary  set  of  the  set  (3.54).  Since  the  set  of  scalings 
used  for  Hoc  gain  scheduling  is  larger  than  that  for  robust  LTI  control,  it  will  typically 
provide  controllers  with  better  performance. 

3.6  Numerical  Issues  of  Controller  Computation 
A  discussion  on  some  of  the  numerical  issues  related  to  solving  for  an  Hoo  gain 
scheduled  controller  would  be  helpful,  since  one  of  the.  major  objectives  is  to  provide 
control  engineers  with  practical  methods  on  designing  gain  scheduled  controllers.  As 
mentioned  earlier  powerful  toolboxes  have  been  developed  to  solve  LMI  equations. 
Two  packages,  in  particular,  are  compatible  with  MATLAB.  One  is  LMI-LAB  [32], 
which  is  a  MATLAB  toolbox.  The  other  is  SDP  Solver  [15]  developed  by  Professor 
S.  Boyd  at  Stanford.  This  package  is  extremely  easy  to  use.  SDP  Solver  accepts 
as  input  and  writes  as  output  .mat  files.  Therefore,  it  can  easily  be  used  inside 
MATLAB  m-files. 

Although  there  are  very  powerful  and  efficient  LMI  solver  packages  available, 
our  particular  problem  tends  to  be  ill-conditioned  and  care  must  be  taken  to  avoid 
numerical  problems.  Even  though  one  of  the  advantages  of  the  LMI  formulation  of 
the  Hoo  problem  is  that  it  relaxes  the  constraint  that  D12  and  D21  must  be  full  column 
and  row  rank  respectively,  for  practical  controller  solutions,  these  matrices  must  be 


87 


full  rank  with  reasonable  values.  Otherwise,  controllers  with  large  poles  will  result 
since  it  does  not  believe  there  is  any  measurement  noise.  The  observer  portion  of  the 
controller  acts  as  a  differentiator,  there  by  resulting  in  an  extremely  fast  controller 
which  is  not  practically  implementable. 

As  one  iterates  solving  the  necessary  and  sufficient  conditions  stated  in  Theorem 
3.5  to  find  7  that  is  close  to  optimal,  the  quadruple  {R,  S,  J3,  L3)  tends  to  have  large 
condition  numbers.  In  fact,  solving  Theorem  3.5  for  large  values  of  7  will  result  in  an 
ill-conditioned  quadruple  (R,  S,  J3,  L3)  if  precautions  are  not  taken.  The  reason  is 
because  the  LMI  solver  stops  as  soon  as  all  of  the  LMI  conditions  are  satisfied.  This 
will  occur  right  on  the  boundary  for  at  least  one  of  the  matrices  in  the  quadruple. 
For  example,  say  that 


R  I 
I  S 


>0 


is  the  active  constraint.  Implied  by  this  LMI  constraint  are  two  more  LMI  constraints: 
R>0  and  5  >  0.  Again  assume  that  i?  >  0  is  the  active  constraint,  then  the  LMI 


solver  will  stop  as  soon  as  it  has  found  an  i?  >  0.  This  R  will  be  near  singular  since 
it  will  be  extremely  close  to  the  boundary.  This  near  singular  R  will  result  in  an  A 
that  is  even  more  ill-conditioned  due  to  the  use  of  matrix  dilation  to  compute  it.  This 
ill-conditioned  X  may  result  in  the  control  designer  not  being  able  to  solve  for  Q,  in 
(3.39)  due  to  numerical  problems. 

To  mitigate  this  problem  the  LMI  solver  needs  to  search  for  a  quadruple  that  is 
not  ill-conditioned  and  yet  satisfies  all  of  the  constraints.  This  will  require  more  than 
just  adding  a  lower  bound  to  (3.29)  and  (3.30)  because  the  condition  number  of  a 
matrix  is  computed  as  the  ratio  of  the  maximum  singular  value  over  the  minimum 
singular  value.  Adding  a  lower  bound  will  not  necessarily  effect  this  ratio  since  the 
maximum  singular  value  is  not  also  constrained.  The  method  by  which  to  constrain 
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the  TnayiTmim  singular  value  is  to  use  the  objective  function  or,  in  terms  of  Linear 
Quadratic  control,  the  cost  function.  If  the  objective  function  is 


R 


Fobj{R,  S,  J3,  L3)  =  Trace 


S 


Jz 


(3.56) 


L  H 

then  the  LMI  solver  will  solve  for  a  quadruple  that  minimizes  this  trace.  Since  the 
trace  of  a  matrix  is  the  sum  of  the  eigenvalues  of  that  matrix,  this  will  minimize  the 
TnflxiTimmn  eigenvalue  which  for  positive  definite  matrices  also  minimizes  the  maximum 
singular  value. 

Another  problem  that  arises  when  solving  these  LMI  constraints  is  that  the  solver 
package  will  return  an  error  indicating  that  the  unknown  variables  are  linearly  depen¬ 
dent.  One  method  to  eliminate  this  problem  is  to  begin  with  a  large  guess  for  7  and 
then  continue  reducing  7  until  it  has  reached  its  optimum  value  with  in  a  tolerance  or 
until  7  <  1  assuming  that  the  parameters  were  scaled  to  have  norm  bound  less  than 
1.  On  the  first  iteration  with  the  large  7  the  initial  guess  for  is  not  crucial  but  on 
successive  iterations  use  the  previous  Q.  for  the  initial  guess.  In  addition,  one  or  more 
elements  of  Q,  may  have  to  be  constrained.  This  can  be  accomplished  by  imposing 
equality  constraints  on  several  elements  of  Q.  The  idea  is  to  constrain  as  few  elements 
as  pdssible  in  order  to  find  a  quadruple.  As  7  is  significantly' reduced  and  good  initial 
guesses  for  Q  are  obtained,  the  number  of  elements  of  that  are  constrained  can 
generally  be  reduced.  When  constraining  elements  of  Cl,  constrain  them  to  have  the 
value  of  the  corresponding  elements  of  Cl  from  the  previous  successful  iteration.  Since 
it  is  desired  to  have  all  of  the  elements  of  Cl  to  be  free  and  not  to  be  constrained,  it 
is  necessary  to  vary  which  elements  are  constrained.  This  has  the  affect  of  allowing 
all  of  the  variables  to  vary  over  several  iterations  but  in  any  one  particular  iteration 
several  elements  may  be  constrained. 


CHAPTER  4 
N-D  Realization 

The  problem  of  realizing  a  linear  parameter-varying  (LPV)  system,  whose  state 
space  model  is  a  rational  function  of  the  elements,  as  an  LFT  is  directly  connected  to 
N-D  system  realization.  N-D  system  realization  has  proven  to  be  a  difficult  area  of 
research  especially  with  respect  to  finding  the  minimal  realization.  In  this  chapter, 
we  begin  by  illustrating  how  realization  theory  may  be  used  to  derive  a  representation 
for  a  LPV  system. 

The  LPV  plant  description  in  terms  of  state  variables  is  given  (see  chapter  2)  as 


x{t) 

m 

B(e) 

x{t) 

yit) 

_c{e) 

0 

u{t) 

In  the  following  we  will  consider  u{t)  as  a  general  input  which  includes  the  dis¬ 
turbance  inputs  (w  in  chapter  2).  We  will  also  consider  that  y{t)  includes  the  usual 
y{t)  outputs  and  the  controlled  or  performance  outputs  q{t). 

Assuming  9  is  constant,  we  may  derive  a  transfer  function  representation  from  (1) 
as 

r(s)  =C(0)(s/-A(0))^^B(0)17(s)  (4.2) 

where  the  transfer  function 


Ci0){sl  -  Aie))-^B{0) 


(4.3) 
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is  a  rational  function  of  s.  Extending  this  idea  to  the  functional  dependence  of  the 
matrices  A{9),  B{0),  and  C{6)  on  9  we  consider  the  relation 

q  =  F{9)p  (4.4) 

and  ask  when  we  can  write 

F{9)  =  Cf{Q~^  -  +  Df  (4.5) 

where  now  we  take  0  as  a  matrix  whose  elements  are  linear  functions  of  the  param¬ 
eters  9.  By  a  slight  abuse  of  notation,  henceforth  9  (without  the  boldface)  will  be 
used  to  refer  to  the  matrix  (previously  0  )  of  parameters.  If  F{9)  can  be  written  as 
(4.5)  exactly  or  approximately,  then  (4.4)  can  be  rewritten  as 

q  =  Cf{9-'^  -  Af)~'^BfP  +  DfP  (4.6) 

Define 

r  =  {9-^-Af)-^BfP  (4.7) 

which  may  be  written  as 

9~^r  =  Aft  +  BfP  (4.8) 

Further  define  w  by  the  relation  r  —  9w  which  leads  to  the  set  of  equations 

q  =  Cft  -I-  DfP  (4.9) 

w  =  Aft  -I-  BfP  (4-10) 

r  —  9w  (4-11) 

Hence  (4.9),  (4.10),  and  (4.11),  which  are  linear  in  9,  may  be  used  to  replace  (4.4) 
with  the  assumption  that  F(9)  is  a  rational  function  of  9. 
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Applying  this  idea  to  the  matrices  in  (4.1)  is  simplified  by  the  definition 


Then 


becomes 


m 


A{d)  Bid) 
0(6)  0 


i(t) 

yit) 


m 


x(t) 

u(t) 


xit) 

yit) 


w 

r 


CpT  +  Dp 


x(t) 

u(t) 


ApT  +  Bp 


x(t) 


dw 


(4.12) 


(4.13) 


(4.14) 

(4.15) 

(4.16) 


In  general  the  matrix  6  may  be  a  linear  function  of  many  parameters.  Therefore, 
assuming  (4.5)  holds,  we  seek  matrices  (Ap,  Bp,  Cp,  Dp)  given  the  matrix  function 
F(6).  Comparison  of  (4.5)  with  the  usual  formula  for  a  transfer  function  of  a  L.T.I 
system  with  matrices  {A,B,C,D) 


C{sl  -  A)-^B  +  D 


reveals  that  si  has  been  replaced  by  d~^.  There  are  two  problems  to  consider  here. 
The  first  is  to  extend  realizations  to  the  case  where  6~^  is  not  a  constant  times  the 
identity.  The  second  is  to  consider  the  selection  of  “small”  matrices  which  work,  i.e.  a 
minimal  realization,  primarily  due  to  the  large  increa.se  in  computational  requirements 
with  the  dimension  of  9.  In  the  following  sections  this  chapter  discusses  non-minimal 
realizations  and  then  turns  to  minimal-realizations. 
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4.1  Non-Minimal  Realizations 

One  possible  realization  of  (4.4)  is  to  realize  each  element  of  the  matrix  F{6) 
separately.  To  illustrate  the  idea,  we  begin  with  a  few  simple  examples. 

To  realize  the  scalar  function 

F(a, h^c)  =  ah  +  c?  -{■  c  (4.17) 

by  inspection  we  construct  the  following  linear  equations; 

w  =  ap  (4.18) 

q  =  bw  +  aw  +  cp  (4-19) 

Then  the  input  p  and  output  q  in  (4.18)  and  (4.19)  realize  the  relationship  in  (4.17). 
For  a  more  involved  example,  suppose  we  wish  to  realize  the  multinomial  matrix 


ab  +  a^  +  c  ac  +  b  x 

ab  +  c  }p  +  a  y 


(4.20) 


The  problem  of  finding  a  minimal  realization  here  is  equivalent  to  computing  p 
and  q  from  x  and  y  with  the  minimum  number  of  multiplications.  As  a  start,  we 
nfight  search  for  a  common  factor  in  the  rows  and  then  search  for  a  common  factor 
in  the  columns.  We  note  there  are  none.  However,  we  can  search  for  common  factors 
in  groupings  of  terms.  The  term 


{ab  +  c)x 
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appears  in  computing  p  and  q  so  it  needs  to  be  computed  only  once.  Furthermore, 
the  term  ax  appears  twice.  Hence  terms  in  the  first  column  may  be  computed  with 


a  =  ax  X 

(4.21) 

Q-first  =  axh  +  cx  X 

(4.22) 

p-first  =  q  +  ax  a 

(4.23) 

where  the  temporary  variable  oc  has  been  introduced.  This  representation  has  four 
multiplications  in  place  of  the  nine  multiplications  which  would  result  with  the  origi¬ 
nal  equations.  The  second  column  terms  may  be  found  and  added  to  the  other  terms 
by 


P 

=  bxy 

(4.24) 

7 

=  axy 

(4.25) 

P 

=  p-first  -H  c  X  7  -f  /? 

(4.26) 

Q 

=  q-first  -t-  6  X  +  7 

(4.27) 

with  the  new  temporary  variables  ^  and  7.  Viewed  in  this  light,  the  reduction 
of  calculations  involved  in  the  non-minimal  realizations  is  related  to  a  problem  in 
optimizing  compilers.  An  optimizing  compiler  attempts  to  find  a  minimal  sequence 
of  operations  to  evaluate  a  given  expression  by  recognizing  and  exploiting  common 
terms  and  by  the  clever  introduction  of  temporary  variables.  It  might  prove  fruitful 
to  see  if  the  techniques  used  by  optimizing  compilers  could  lead  to  algorithms  for 
minimal  realizations. 

It  is  clear  to  see  how  we  can  realize  any  F(9)  which  is  a  multinomial  in  9.  If  F{6) 
is  a  scalar  rational  function  of  9, 

n(») 
die) 


Fie) 


(4.28) 
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where  n(6)  is  a  multinomial  in  9  and  d{9)  is  of  the  form  1  —  D{9)  where  D{6)  is  a 
multinomial  in  9,  then  a  simple  technique  works. 

We  may  realize  n{9)  and  D{9)  =  1  -  d{9)  as  before  and  then  construct  F{9)  by 
the  feedback  formula  as  follows: 

q  =  e  (4.29) 

e  =  n{9)p-Di9)q  (4.30) 

Currently  minimal  realizations  exist  for  the  following  special  cases:  a)  2-D  transfer 
functions  with  separable  numerator  [37],  [13];  b)  2-D  transfer  functions  with  separable 
denominator  [38],  [13];  c)  2-D  all-pole  and  all-zero  [39],  d)  2-D,  3-D,  and  N-D  systems 
that  can  be  expanded  into  a  continued  fraction  expansion  (CFE)  [40],  [41],  [42],  [43], 

[44],  [45].  It  should  be  noted  that  all  of  these  special  case  methods  apply  to  the  single¬ 

input  single-output  case  (SISO)  systems.  For  a  missile  autopilot  design  and  for  many 
other  practical  systems  a  multi-input  multi-output  (MIMO)  method  is  required. 

In  [46]  Koonar  and  Sohal  present  a  more  general  method  for  2-D  realization  us¬ 
ing  Markov  parameters  and  a  moments  of  impulse  response  matrix.  In  [47]  Eydgahi 
provides  a  counter  example  to  Koonar  and  Sohal’s  method.  After  searching  the  liter¬ 
ature,  no  response  has  been  discovered  to  this  counter  example.  Another  limitation  to 
Koonar  and  Sohal’s  method  is  that  it  requires  that  the  total  degree  of  the  numerator 
polynomial  should  at  most  be  equal  to  that  of  the  denominator  polynomial.  This  is 
a  drawback  for  our  case  since,  in  many  instances,  linear  multi-parameter-systems  are 
curve  fitted  with  N-D  polynomials.  This  results  in  LPV  systems  whose  elements  are 
multivariate  polynomials  in  the  parameters,  for  which  Koonar  and  Sohal’s  method 
would  not  be  valid.  Recently,  general  nonminimal  realization  methods  have  been 
developed  for  both  SISO  and  MIMO  N-D  systems. 
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4.2  General  Nonminimal  Realization  Methods 
Lambrecht’s  paper  [48]  presents  a  general  method  by  which  to  compute  an  LFT 
system  description  of  a  linear  system  with  parametric  uncertainties.  The  entries  of  the 
state-space  matrices  are  bounded  and  can  be  given  as  real-rational  N-D  polynomials 
in  the  parameters.  His  method  begins  by  modeling  the  most  basic  element  as  an 
LFT,  i.e.  an  individual  term  of  an  N-D  polynomial.  Once  the  most  basic  element  is 
modeled  they  then  exploit  the  fact  that  parallel  combinations  of  LFTs,  cascading  of 
LFTs,  and  inverses  of  LFTs,  all  result  in  LFTs.  From  these  properties,  they  develop 
one  LFT  for  the  entire  system.  This  LFT  will  be  quite  large  since  this  method  does  not 
take  into  account  any  redundancy.  They  then  reduce  this  large  LFT  by  treating  the 
parameters  individually  as  states  and  then  extracting  out  the  locally  “uncontrollable” 
and  “unobservable”  parameters.  This  will  reduce  the  LFT  substantially  but  it  is  not 
clear  how  close  the  resulting  LFT  will  be  to  minimal. 

Cheng  and  DeMoor  [49],  [50]  provide  an  algorithm  to  systematically  solve  the  N- 
D  realization  problem.  They  do  not  provide  a  minimal  realization  but  a  systematic 
method  to  obtain  a  realization.  The  authors  do  suggest  how  to  reduce  the  order  of 
the  realization  which  requires  the  solution  of  a  set  of  matrix  equations  such  that  the 
solution  matrix  has  minimal  rank.  Unfortunately,  it  is  still  an  open  question  on  how 
to  find  the  minimal  rank  solution  to  this  set  of  equations. 

Cheng  and  DeMoor  consider  the  following  N-D  realization  problem.  Given  a 
multivariate  rational  matrix  F{9)  E  72-"^"*,  find  an  LFT  with  block  structure  bg  = 
[ri,r2, . . .  jTg]  with  being  an  integer  and  a  coefficient  matrix  L  E  7:2,(”+Ox ("»+»■) ^ 
r  =  Yli=i  partitioned  as 


(4.31) 
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such  that 


F{9)  =  L22  +  L2i9{I  -  9Ln)~^Li2  (4.32) 

Cheng  and  DeMoor  assumed  the  structure  for  9  was  diagonal 

9  =  Siag{9Jr^^■■0gIr,)  (4.33) 

The  elements  of  F{9)  are  further  assumed  to  be  of  the  following  form 

[f  (»)!«  =  ci,  +  Me)  (4.34) 

with  Cij  being  constant  and  fij{9)  rational  multinomials  of  the  form 

^.M0)=<(0)=0.  (4.33) 

and  k(9)  and  l{9)  are  multinomials  in  the  elements  of  9.  Given  the  form  (4.34),  F{9) 
can  be  written  as 

I  '  ■  ■  . 

F(^9)  =  Fo  +  (4.36) 

i— 1 

where  Fq  is  a  constant  matrix  with  entries  Cij.  The  functions  fi{9),  i  =  1, . .  .m  are 
formed  from  the  elements  fij{9)  of  F{9)  such  that  for  each  pair  i,j  there  is  a  A:  such 
that  fk(9)  =  fijid).  Finally  Fi  is  the  constant  coefficient  matrix  that  tells  how  fi{9) 
enters  F{9). 

The  algorithm  for  constructing  a  realization  starts  by  first  realizing  (4.36)  which 
is  a  linear  function  of  the  /i(0)’s.  That  is  we  find  a  coefficient  matrix  La  so 


F{9)  —  L022  +  La2i^{f){I  -  ^Lan 


(4.37) 
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where 

La  =  ,  K  =  ,ra2,---,  %]  (4.38) 

■^a21  ^(122 

A(/)  =  diag(/,(«)4., . (4.39) 

Safanov  and  Athans’  Internal  Feedback  Loop  ([51])  parameter  representation  is 
used  to  find  La-  A  realization  with  =  0  can  be  easily  found  by  first  factoring 
each  of  the  Fi  matrices  in  (4.36)  as 

Fi  =  Ui'LiVT 

by  singular  value  decomposition  with  Sj  a  Voi  x  ra.  positive  definite  diagonal  matrix. 
Then  if 

U 

F(0) 

Comparison  with  (4.37)  gives 

La22  ~  ^0 
La2^  =  U 
La,2  = 

La,,  =  0  (4.41) 

The  next  step  involves  finding  a  realization  for  A(/)  as 


=  [E?T^^,Etvf,...,EfV[^] 

=  Fo  +  t/A(/)V^  (4.40) 


A(/)  =  ite  +  it,,  A(«)(/  -  it..  A(«))-‘Lt., 
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Assume  that  each  (scalar)  rational  function  fi{9)  has  a  realization  Then  the 
realization  Lb  above  can  be  obtained  by  combining  these  realizations.  The  bookkeep¬ 
ing  for  this  procedure  is  a  bit  tedious  but  possible.  By  the  assumption  /i(0)  =  0, 
A{f)ff=o  =  0  so  that  Lb22  =  0.  The  final  answer  becomes 

Afcii  -^612  ■^012 

■^021  ■^621  ^0,22 

It  is  clear  that  this  realization  is  by  no  means  minimal,  however  Cheng  and  De- 
Moor  give  some  suggestions  for  reducing  the  complexity. 

In  a  recent  paper,  Belcastro  et.  al.  [11]  present  an  overview  of  a  matrix-based 
method  for  low-order  realization  of  a  plant  that  has  a  multivariate  rational  functional 
dependence  on  uncertain  parameters.  The  full  details  of  this  method  are  given  in  Bel- 
castro’s  dissertation  [52].  The  claim  of  the  authors  is  that  this  method  is  efficient  and 
produces  a  lower  order  realization  than  other  methods  that  use  reduction  techniques 
such  as  removing  unobservable  and  uncontrollable  modes. 

Belcastro,  et.  al.,  used  the  same  assumptions  and  problem  formulation  presented 
earlier  with  the  exception  that  the  F(0)  matrix  was  given  by 


F{e)  = 


A(0)  B(0) 

cm  Dm 


and  the  parameter  vector  0  was  defined  in  terms  of  a  nominal  parameter  value  as 


9  —  ^0  +  ^ 


Although  Belcastro,  et.al.  do  not  propose  to  find  a  minimal  solution,  they  attempt 
to  find  a  realization  with  a  small  dimension.  Expanding  the  assumed  form  for  F(9) 
in  terms  of  the  nominal  parameters  9o  and  a  perturbation  S  3delds  the  following  LFT 
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formulation  for  F{9). 


F{e)  =  P22  +  P2i{I  -  AP12  =  FiOo)  +  F^iS)  (4.42) 


where 


Fa{S)  =  L2i{I  -  AiS)Ln)-^Ai5)Lu  (4.43) 


If  the  state  matrices  A{0),  B{9),  C{9)  and  D{9)  are  multinomial  functions  of  9. 
straightforward  algebra  can  be  used  to  determine  the  matrices  in  Fa{S).  Prom  (4.42) 
we  have 


m 


A(e)  B{e) 

c())  D(e)_ 

A(9o  +  5)  B{9o  +  5) 
C{9o  +  5)  D{9o  +  5) 


(4.44) 

(4.45) 


which  may  be  expanded  to  separate  the  nominal  and  perturbed  parts  of  the  system. 

In  order  to  equate  (4.43)  with  the  perturbed  part  of  F{9)  from  (4.45),  Fa{S)  must 
be  a  multinomial  function  of  9  which  implies  that  the  expension  of  Fa  (S)  terminates 
or  A((5)Lii  is  nilpotent  with  index  of  nilpotency  r, 


(A(S)lnr*'  =  0 


and 


(I  —  ALii)  ^  =  /  +  (ALii)  +  (ALii)^  +  . . .  (ALii)’'  (4.46) 

However  if  the  state  matrices  are  rational  multinomial  functions  the  expansion 
becomes  more  difficult.  In  her  dissertation  [52],  Belcastro  proposed  first  representing 
the  system  in  a  matrix  fraction  description  MFD  and  presents  a  scheme  for  the 
factorization  of  a  rational  Sa(<S)  as  a  matrix  fraction. 
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SM  =  (4.47) 

where  Sna(^)^  'S'da(<^)  ^d^S)  and  SnA^)  multinomial  matrices. 

By  substituting  (4.46)  into  (4.43),  and  then  comparing  the  result  with  the  Sa{S) 
computed  from  the  MFD  and  then  broken  out  such  that  like  terms  can  be  com¬ 
puted,  results  in  a  set  of  equations  that  can  then  be  used  to  solve  for  Ln,  L12  and 
L21.  If  a  solution  does  not  exist  then  the  order  of  the  offending  5i  is  increased  by 
1  until  a  solution  is  found.  By  increasing  the  order  of  A  only  when  it  is  necessary 
to  obtain  a  solution,  ensures  that  a  low-order  realization  is  found.  It  is  not  clear, 
however,  whether  or  not  the  solution  takes  into  account  any  redundancy  or  possible 
cancellations  between  and  Sd^. 

4.3  Minimal/Near  Minimal  2-D  Realization 
Currently  two  methods  of  realization  show  promise  in  providing  minimal/near 
minimal  2-D  realization.  Realization  methods  for  1-D  systems  have,  as  an  underly¬ 
ing  principle,  the  fact  that  a  minimal  realization  is  both  observable  and  controllable. 
Therefore  one  can  take  a  nominal  1-D  system,  test  for  controllability  and  observabil¬ 
ity  and  then  obtain  a  minimal  realization  by  throwing  away  the  uncontrollable  and 
unobservable  states.  Kung  et.al.  [13]  showed  that  the  notions  of  controllability  and 
observability,  although  extending  to  2-D  systems,  do  not  provide  a  test  for  minimal¬ 
ity.  In  fact  Kung  et.al.  provide  an  example  of  a  controllable  and  observable  system 
that  is  not  minimal.  In  their  paper,  they  investigate  several  ideas  on  determining  if 
a  system  is  minimal.  In  the  end  they  extend  Rosenbrock’s  minimality  test  to  2-D 
systems  which  will  also  hold  for  N-D  systems  as  well.  This  fesult  states  that  a  plant 


p  =  ce[i-Ae]-^B 
=  C[e-^-A]~^B, 


(4.48) 
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where  0  is  as  defined  previously,  is  minimal  if  and  only  if  [0“^  —  A] ,  B  are  left 
coprime  and  C,  [0“^  -  A]  are  right  co-prime. 

There  are  two  difficulties  with  this  definition  of  minimality.  One  is  that  it  allows 
the  A,B  and  C  matrices  to  have  complex  gains  and  the  other  is  that  Youla  and 
Gnavi  [53]  state  that  certain  decomposition  techniques,  which  have  proven  basic  for 
1-D  and  2-D  systems,  are  no  longer  applicable  for  3-D  and  larger  systems.  Therefore 
finding  techniques  to  compute  greatest  common  divisors  (GCD)  for  3-D  and  larger 
systems  will  be  diflficult.  For  2-D  systems  though  Morf  et.  al.  [12]  have  developed 
methods  to  compute  greatest  common  right  (left)  divisors  (GCRD/GCLD).  We  have 
implemented  Morf  et.al’s  GGD  extraction  method  in  MAPLE®.  The  difficulty  that 
Kung  et.al.  had  in  implementing  the  extraction  of  GCDs  to  reduce  non-minimal  state 
descriptions  to  minimal  ones,  was  that  after  the  GCD  had  been  extracted  the  plant 
description  was  no  longer  in  “state  space”  form  (i.e.  parameters  only  on  the  diagonal) 
and  they  did  not  know  how  to  return  it  to  “state  space”  form. 

We  discovered  that  Rosenbrock  in  [14]  introduced  the  concepts  of  the  system 
matrix  and  system  equivalence  for  1-D  systems.  Rosenbrock  defined  for  1-D  systems 
a  system  matrix  by  writing  the  Laplace  transformed  state  equations  in  the  following 
form. 

T{s)  U{s) 

-l/(s)  W(s) 

' - v - 

5 

where  S  is  defined  as  the  system  matrix.  Note  that  the  dimension  of  r(s)  is  n  x  n 
Two  system  matrices  are  system  equivalent  if  and  only  if  they  give  rise  to  the  same 
transfer  function  matrix.  The  relation  of  system  equivalence  is  generated  by  the 
following  operations  on  system  matrices. 

(i)  Multiply  any  one  of  the  first  n  rows(columns)  by  a  rational  function  not  iden- 


-a;(s) 

0 

u{s) 

y{s) 

(4.49) 


tically  zero. 
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(ii)  Add  a  multiple,  by  a  rational  function,  of  any  one  of  the  first  n  rows(columns) 
to  any  other  row  (column). 

(iii)  Interchange  any  two  among  the  first  n  rows  (columns). 

1  0 

(iv)  Add  a  row  and  a  column  to  S  to  form 

0  S 

System  equivalent  operations  can  be  used  to  return  the  system  matrix  back  into  state- 
space  form  after  the  GCD  has  been  extracted.  These  concepts  and  operations  extend 
readily  to  N-D  systems  by  replacing  s  with  0“^.  The  purpose  for  replacing  s  with 
0“^  has  to  do  with  the  nature  of  LFTs.  For  a  linear  system,  the  integrators  (s“^) 
enter  in  a  linear  fractional  manner  and  for  our  LPV  system,  0  and  not  0“^  enters 
in  a  linear  fractional  manner.  Therefore  s  and  0“^  are  similar.  Compare  (4.48)  with 
C[sI-AY^B. 


By  performing  system  equivalent  operations,  one  can  put  the  system  matrix  back 
into  state-space  form.  Below  is  an  algorithm  to  obtain  a  2-D  minimal/near  minimal 
realization. 

Algorithm  4.1  (GCD  Based).  . 

1.  Obtain  a  nonminimal  description  of  the  system. 

2.  Using  [12]  2-D  GCD  extraction  method,  extract  the  GOLD  of  [0~^  —  A],  B. 
If  these  are  already  coprime  then  extract  the  GCRD  of  C,  [0“^  —  A\.  After 
extraction  the  system  matrix  is  no  longer  in  “state  space”  form.  Its  new  general 

T  U 

form  is  notated  as  S  = 

-V  W 
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3.  Use  system  equivalent  operations  to  return  the  T  matrix  back  into  “state  space” 
form  (i.e.  parameters  only  on  the  diagonal).  This  mil  yield  a  new  system 


matrix  to  be  S  = 


Ti  C/i 
-Vi  Wi 


4.  Solve  VgTi^Us  —  ViTf^Ui  —  =  0  for  Vs  and  Ug  such  that  they  contain  only 

real  gains.  Note  that  Tx  is  already  in  “state  space”  form  while  Ux,V\  and  Wx 
are  most  likely  not. 


5.  If  no  real  gain  solution  exists  then  introduce  another  parameter  and  expand  the 
Tx  matrix  up  one  dimension. 


6.  Repeat  steps  4  and  5  until  a  real  gain  solution  is  obtained.  The  real  gain  solution 
will  be  the  minimal/near  minimal  solution  that  is  desired. 


A  difficulty  is  that  the  GCD  extraction  process,  in  general,  reduces  the  system 
too  far  since  it  allows  for  complex  gains.  For  example,  there  are  2-D  systems  that  are 
real  gain  minimal  but,  when  Rosenbrock’s  minimality  test  is  applied,  are  not  found 
to  be  minimal.  The  reason  is  that  there  exist  a  complex  gained  system  that  is  of 
smaller  order  than  the  real  gain  system.  This  is  a  serious  problem  in  the  search  for 
minimal  realization  algorithms  for  N-D  systems.  The  above  algorithm  manages  this 
problem  in  step  5  by  expanding  up  another  dimension.  The  key  here  is  on  intelligently 
expanding  so  that  you  can  get  step  4  to  solve  with  a  real  gain  solution  with  as  few 
of  parameters  as  possible.  Work  is  currently  being  conducted  on  how  to  best  expand 
the  system. 
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4.3.2  System  Equivalent  Based  Algorithm 

The  second  method,  which  shows  promise  for  the  general  case,  is  to  begin  with 
the  plant  P  in  the  lower  right  corner  of  Rosenbrock’s  system  description  matrix 

I  0 
0  P 

The  idea  is  to  move  P  over  to  the  left  using  system  equivalent  operations. 


I  I 
-P  0 

Prom  here  we  can  multiply  the  column  corresponding  to  the  last  column  of  P  by 
the  least  common  multiple  of  the  denominators  of  P  and  then  expand  it  up  until  all 
parameters  are  only  on  the  diagonal  of  T. 

T  U 
-V  W 

Once  T  only  has  parameters  on  the  diagonal  it  is  time  to  begin  eliminating  the 
parameters  out  of  U  and  V.  At  this  point  we  can  use  step  4  of  the  GCD  based 
algorithm  to  see  if  there  exists  a  real  gain  solution.  If  there  is  not  then  we  continue 
expanding  until  a  real  gain  solution  exists.  The  following  is  an  algorithm  for  2-D 
systems  that  implements  this  method. 

Algorithm  4.2  (System  Equivalent  Based).  . 

where  z  and  w  are  the  parameters.  H{z,  w) 

must  have  a  constant  term  in  its  collective  denominator  if  it  does  not  then  create 
one  by  a  change  of  variable. 

2.  Substitute  z  =  —  and  w  —  —  into  Hiz,  w). 

Zi  Wi  \  / 


1.  LetH{z,w)  = 


A{z,w)  B{z,w) 
C{z,w)  D{z,w) 
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3.  Place  H{zi,Wi)  into  the  W  position  of  Rosenbrocks  system  matrix. 

4.  By  use  of  system  equivalent  operations  move  H{zi,Wi)  over  one  position  to  the 
left. 

5.  Clear  the  fraction  of  the  last  column  of  H{zi,  Wi)  by  multiplying  the  correspond¬ 
ing  column  by  d{zi,  Wi)  which  is  the  least  common  multiple  of  denominators  of 
H{zi,Wi).  Now  the  T  matrix  contains  the  least  common  multiple  of  the  denom¬ 
inators. 

6.  Expand  the  T  matrix  up  until  it  is  in  “state  space”  form.  This  can  be  done  by 
moving  d{zi,  Wi)  over  to  the  left  one  space  and  at  the  same  time  depositing  a  Zi 
in  its  place  on  diagonal.  Then  by  column  multiplication  eliminate  all  of  the  Zi ’s 
from  the  d(zi,Wi).  This  will  place  a  new  di{zi,Wi)  with  the  order  of  Zi  reduced 
by  one  on  the  diagonal  above.  Continue  depositing  Zi ’s  until  they  are  completely 
eliminated  and  then  by  the  same  process  place  Wi ’s  on  the  diagonal  until  they 
are  eliminated. 

7.  Solve  VsT~Ws  —  VT~^U  —  W  =  0  for  Vg  and  Ug  such  that  they  contain  only 
real  gains.  Note  that  T  is  already  in  “state  space”  form  while  U ,  V  and  W  are 
most  likely  not. 

8.  If  no  real  gain  solution  exists  then  introduce  another  parameter  and  expand  the 
Ti  matrix  up  one  dimension. 

9.  Repeat  steps  7  and  8  until  a  real  gain  solution  is  obtained.  The  real  gain  solution 
will  be  the  minimal/near  minimal  solution  that  is  desired. 

Due  to  the  special  structure  that  this  algorithm  creates  by  placing  I’s  on  the 
super-diagonal,  it  is  not  completely  general  for  finding  the  minimal  realization  but 
can  find  near  minimal  realizations  for  all  systems. 
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The  idea  of  starting  with 


matrix  in  the  lower  right  of  the 


(7(0-1)  £>(0-1) 

system  matrix  and  then  expanding  up  can  easily  be  extended  to  N-D  systems.  The 
above  algorithm  will  always  put  T  in  “state  space”  form  for  2-D  systems  but  for  N-D 
systems,  in  general,  it  will  not. 

Both  of  these  methods  work  on  both  MIMO  and  SISO  systems.  The  second 
method  shows  promise  on  obtaining  near  minimal  realization  for  general  N-D  systems. 


CHAPTER  5 
Conclusion 

In  this  report,  we  focused  on  the  problem  of  gain  scheduling  LPV  systems.  In  par¬ 
ticular  we  developed  necessary  and  sufficient  conditions  for  determining  the  Quadratic 
stabilizability  of  LPV  systems  and  provide  a  stabilizing  LPV  controller  if  one  exists. 
Our  necessary  and  sufficient  conditions  are  realistically  computable  with  reduced  di¬ 
mension  and  with  only  the  elements  of  the  Lyapunov  P  as  the  unknown  variables 
where  as  the  one  given  by  Becker  in  [25]  is  of  full  dimension  and  requires  an  addi¬ 
tional  set  of  variables  the  size  of  the  controller  for  each  grid  point.  Examples  were 
given  that  demonstrated  problems  that  can  be  encountered  with  “frozen  coefficient” 
gain  scheduling  which  can  be  overcome  by  our  method.  We  investigated  the  geomet¬ 
ric  interpretation  of  the  A22  condition  and,  for  the  second  order  case,  developed  a 
method  by  which  to  graphically  display  the  set  of  all  stabilizing  Lyapunov  functions 
for  an  LPV  system.  Furthermore  our  solution  to  the  A22  condition  is  usefull  in  solving 
Becker  and  Packards  [6]  output  stabilization  problem  for  which  they  set  up  but  did 
not  provide  a  solution. 

We  extended  our  results  on  Quadratic  stability  to  the  more  general  case  of  General 
Lyapunov  stability.  The  derivation  was  the  same  as  for  Quadratic  stability,  except 
there  was  an  additional  term  that  was  a  function  of  P.  We  developed  a  method  by 
which  to  solve  these  new  conditons.  In  fact,  when  Quadratic  stabilizability  fails,  we 
showed  that  it  was  necessary  to  split  the  parameter  space  into  regions  of  frozen  Lya¬ 
punov  functions.  After  determining  all  of  the  necessary  frozen  Lyapunov  functions, 
these  functions  are  then  “stitched”  together  such  that  General  Lyapunov  stability 
holds.  This  requires  knowledge  of  the  rate  of  change  of  the  parameter  variations. 
In  order  to  properly  stitch  these  Lyapunov  functions  together,  P  cannot  change  so 
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fast  that  the  necessary  and  sufficient  conditions  for  General  Lyapunov  stability  no 
longer  hold.  This  fact  is  taken  into  account  when  first  determining  the  frozen  regions. 
Now,  instead  of  the  clasical  frozen  point  controllers,  we  have  frozen  region  Lyapunov 
functions.  Prom  these  Lyapunov  functions  we  showed  how  to  construct  a  stabilizing 
controller. 

In  addition,  we  add  performance  to  our  Quadratic  stability  results  for  strictly 
proper  systems  by  guaranteeing  that  the  closed  loop  system  has  an  induced  C2  norm 
less  than  a  specified  value  if  a  modified  A22  condition  holds.  A  missile  autopilot 
example  was  provided. 

Another  method  of  gain  scheduling  LPV  systems  is  via  Hoc  based  gains  scheduling 
as  developed  by  Apkarian  and  Cabinet  [8].  We  reviewed  their  work  and  discussed 
how  to  construct  a  controller  and  how  to  mitigate  the  inherent  numerical  problems 
associated  in  solving  the  LMFs.  This  method  of  LPV  stabilization  does  not  require 
gridding  the  parameter  space  but  also  does  not  take  into  account  the  realness  of  the 
parameters  as  does  our  Lyapunov  based  method  of  gain  scheduling.  It  also  requires 
that  a  realization  of  the  plant  with  a  linear  fractional  dependence  on  the  parameters 
be  provided. 

In  the  final  chapter  we  discussed  methods  of  finding  linear  fractional  realizations. 
These  realizations  are  related  to  N-D  realizations  for  which  minimal  realizations  are 
very  difficult  to  find.  We  developed  two  methods  for  finding  minimal/near  minimal 
realizations  for  2-D  systems  based  on  the  work  of  Kung  et.al.  [12, 13]  and  Rosenbrock 
[14].  The  first  method  required  extracting  any  greatest  common  divisors  out  of  a 
nonminimal  realization  and  then  returning  the  realization  back  to  "state-space”  form. 
This  method  does  not  extend  to  N-D  systems  since  the  techniques  used  to  compute 
the  GCD  for  2-D  sytems  does  not  extend  to  N-D  systems  [53].  However  the  second 
method,  the  System  Equivalent  Based  method,  can  be  extended  to  N-D  systems.  This 
method  begins  with  the  LPV  system  description  being  converted  to  Rosenbrock’s 
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system  matrix  form.  Then  the  system  matrix  form  is  systematically  manipulated 
until  it  is  in  “state-space”  form. 

This  report  provides  the  control  engineer  with  a  design  method  of  gain  scheduling 
LPV  systems  that  relaxes  the  requirement  of  scheduling  on  slowly  varying  parame¬ 
ters.  Missile  systems  have  been  approximated  as  LPV  systems  in  [54]  [55]  and  [28] 
to  list  a  few  examples.  We  demonstrated  in  Section  2.3  that  this  work  can  now  be 
applied  to  the  missile  autopilot  problem.  In  addition,  there  are  many  other  succes- 
ful  applications  of  classical  gain  scheduling  such  as  aircraft  control  systems,  turbojet 
engines,  ship  steering,  pH  control,  combustion  control  and  fuel-air  control  in  an  au¬ 
tomobile  engine  (see  for  example,  [56,  57,  58])  for  which  our  method  of  Lyapunov 
based  gains  scheduling  could  be  applied. 
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